Preprint typeset using lAT^^X style cmulateapj v. 11/12/01 



AN EXPANDED RXTE SURVEY OF X-RAY VARIABILITY IN SEYFERT 1 GALAXIES 



A. Markowitz"''^, R. Edelson 



.t3 



^ Laboratory for High-Energy Astrophysics, NASA/Goddard Space Flight Center, 
Code 662, Greenbelt, MD 20771; agm@milkyway.gsfc.nasa.gov 
^ N.A.S./N.R.C. Research Associate 
^ Dept. of Astronomy, Univ. of California, Los Angeles CA 90095-1562 



o 
o 

(N 

CD 



> 

in 

o 

o\ 
o 

o 

:^ 

Oh' 

6 



X 



ABSTRACT 

The first seven years of RXTE monitoring of Seyfert 1 active galactic nuclei have been systematically 
analyzed to yield five homogeneous samples of 2-12 keV light curves, probing hard X-ray variability 
on successively longer durations from ~1 day to '^3.5 years. 2-10 keV variability on time scales of ~1 
day, as probed by ASCA, are included. All sources exhibit stronger X-ray variability towards longer 
time scales, but the increase is greater for relatively higher luminosity sources. Variability amplitudes 
are anti-correlated with X-ray luminosity and black hole mass, but amplitudes saturate and become 
independent of luminosity or black hole mass towards the longest time scales. The data are consistent 
with the models of power spectral density (PSD) movement described in Markowitz et al. (2003a) and 
McHardy et al. (2004), whereby Seyfert 1 galaxies' variability can be described by a single, universal 
PSD shape whose break frequency scales with black hole mass. The best-fitting scaling relations between 
variability time scale, black hole mass and X-ray luminosity imply an average accretion rate of ~5% of 
the Eddington limit for the sample. Nearly all sources exhibit stronger variability in the relatively soft 
2-4 keV band compared to the 7-12 keV band on all time scales. There are indications that relatively 
less luminous or less massive sources exhibit a greater degree of spectral variability for a given increase 
in overall flux. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

X-ray observations can provide constraints on the phys- 
ical conditions in the innermost regions of Seyfert 1 Ac- 
tive Galactic Nuclei (AGNs), as the X-rays are generally 
thought to originate in close proximity to the putative cen- 
tral supermassive black hole. On the basis of spectroscopic 
observations, the leading models of the X-ray continuum 
production are based on a hot, Comptonizing electron or 
electron-positron pair corona close to the black hole. The 
exact geometry remains uncertain, though numerous mod- 
els have been invoked (e.g, Zdziarski et al. 2003), includ- 
ing a neutral accretion disk extending in to the minimum 
stable orbit and sandwiched by a patchy and possible out- 
flowing corona (e.g.. Stern et al. 1995, Svensson 1996, Be- 
loborodov 1999) and a hot inner disk radially surrounded 
by a cold disk, with a variable transition radius (Shapiro, 
Lightman & Eardley 1976, Zdziarski, Lubinski & Smith 
1999). The corona multiply- upscatters thermal soft pho- 
tons emitted from the disk to produce an X-ray power-law 
in the energy range 1-100 keV (e.g., Haardt, Maraschi 
& Ghisellini 1994). Furthermore, the disk, or some other 
cold, optically thick material, reprocesses the hard X-rays, 
as evidenced by the so-called 'Compton reflection humps' 
above ~10 keV in Seyfert spectra, as well as strong iron 
fluorescent lines at ~6.4 keV (Lightman & White 1988, 
Guilbert & Rees 1988, Pounds et al. 1990). 

Seyfert 1 galaxies exhibit rapid, aperiodic X-ray con- 
tinuum variability for which no fully satisfying explana- 
tion has been advanced. Probably the best way to char- 
acterize single-band Seyfert variability, if adequate data 
exist, is to measure the fluctuation power spectral density 



(PSD) function. Recent studies such as Edelson & Nandra 
(1999), Uttley, McHardy & Papadakis (2002), Markowitz 
et al. (2003a), Marshall et al. (2004) and McHardy et 
al. (2004) measured high-dynamic range broadband PSDs 
which showed the red-noise nature of Seyfert variability 
at high frequencies, but flattened below temporal frequen- 
cies corresponding to time scales of a few days. Markowitz 
et al. (2003a) developed a scenario in which all Seyfert Is 
have a PSD shape similar to that of X-ray Binaries (XRBs) 
and which scale towards lower temporal frequency with in- 
creasing black hole mass. Physically, this is consistent with 
a scenario in which relatively more massive black holes 
host larger X-ray emitting regions, the variability mecha- 
nism takes a longer time to propagate through the emission 
region, and the observed variability is 'slower.' 

When data are not adequate to construct a PSD, it is 
still valuable to quantify the variability amplitude. The 
well-known anticorrelation between variability amplitude 
(as quantified over a fixed temporal frequency range) and 
source luminosity on both short time scales ('^1 d: Barr & 
Mushotzky 1986; Nandra et al. 1997a, Turner et al. 1999) 
as well as long time scales ('^300 d: Markowitz & Edel- 
son 2001, hereafter MEOl) is consistent with the above 
physical interpretation. 

Numerous X-ray spectral variability studies (e.g., 
Markowitz, Edelson & Vaughan 2003b; also Nandra et 
al. 1997a, MEOl) have shown the majority of Seyferts to 
soften as they brighten, with the relatively softer ener- 
gies displaying stronger variability. It is currently unclear 
whether this is due to intrinsic slope changes of the coronal 
power-law continuum (e.g.. Lamer et al. 2003a, Uttley et 
al. 2003) or due to the presence of a much less variable hard 
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component that is likely associated with the Compton re- 
flection hump (e.g. Shih, Iwasawa & Fabian 2002; Taylor, 
Uttley & McHardy 2003). In contrast to the 'normal' or 
'broad-line' Seyfert I's which show this property, however, 
some 'Narrow-Line' or 'soft-spectrum' Seyfert Is (charac- 
terized by FWHM < 2000 km s~^, and steep photon in- 
dices; e.g., Boiler, Brandt & Fink 1996) have been seen to 
vary with a much weaker dependence on energy compared 
to broad-line Seyfert Is (e.g., Edelson et al. 2002, Vaughan 
et al. 2002). 

The archival data accumulated by the Rossi X-ray Tim- 
ing Explorer [RXTE) during its first seven years of opera- 
tion permits a study of broadband continuum and spectral 
variability behavior on time scales ranging from days to 
years. The long-term variability survey of MEOl was the 
first to systematically probe X-ray variability on such long 
time scales, examining nine Seyfert 1 light curves each of 
300 days in duration. This paper expands that survey to 
cover additional time scales and sources using additional 
archival RXTE data. In this paper we test the relation be- 
tween X-ray variability and black hole mass, including the 
idea of broadband PSD movement with black hole mass, 
and exploring spectral variability throughout Seyfert Is. 
The source selection and data reduction are described in 
§2. The sampling and analysis are described in §3. The 
results are discussed in §4, and a short summary is given 
in §5. An Appendix briefly explores if the modeled RXTE 
PCA background has any significant effect on the mea- 
sured variability properties for low count rate or steep- 
spectrum sources. 

2. DATA COLLECTION AND REDUCTION 

RXTE has observed ~55 Seyfert 1 galaxies during the 
first seven years of its mission. Data taken through most 
of Cycle 7 had turned public by 2004 February, when these 
analyses were performed. This paper considered these 
data as well as the authors' proprietary observations of 
three Seyfert 1 galaxies observed during Cycle 8. §2.1 de- 
tails how the RXTE data were reduced. 

The observational approach of this project was to obtain 
monitoring on multiple long time scales, sampled as uni- 
formly as possible for as many Seyfert 1 galaxies as possi- 
ble. Using the available archive of RXTE data to optimize 
this trade-off yielded a sample of 27 Seyfert Is suitable for 
analysis on at least one of the time scales of interest, 1 d, 
6 d, 36 d, 216 d, or 1296 d. Additionally, most of these 
sources also had adequate short time scale (1 d) ASCA 
data publicly available. Most of the sources with data on 
the 36 d, 216 d, and 1296 d time scales have had their PSDs 
measured or are currently undergoing monitoring for fu- 
ture PSD measurement. §2.2 and §2.3 detail construction 
of the RXTE and ASCA light curves, respectively. 

2.1. RXTE data reduction 

All of the RXTE data were taken with the Proportional 
Counter Array (PCA), which consists of five identical col- 
limated proportional counter units (PCUs; Swank 1998). 
For simplicity, data were collected only from those PCUs 
which did not suffer from repeated breakdown during on- 
source time (PCUs 0, 1, and 2 prior to 1998 December 23; 
PCUs and 2 from 1998 December 23 until 2000 May 12; 

* see also http://lheawww.gsfc.nasa.gov/users/craigm/pca-bkg/bkg-users 



PCU 2 only after 2000 May 12). Count rates quoted in this 
paper are normalized to 1 PCU. Only PCA STANDARD-2 
data were considered. The data were reduced using stan- 
dard extraction methods and FTOOLS v5.2 software. 
Data were rejected if they were gathered less than 10° 
from the Earth's limb, if they were obtained within 30 min 
after the satellite's passage through the South Atlantic 
Anomaly (SAA), if ELECTRONO > 0.1 (ELECTRON2 
after 2000 May 12), or if the satellite's pointing offset was 
greater than 0?02. 

As the PCA has no simultaneous background monitor- 
ing capability, background data were estimated by using 
PCABACKEST v2.1e to generate model files based on the 
particle-induced background, SAA activity, and the diffuse 
X-ray backgroiuid. This background subtraction is the 
dominant source of systematic error in RXTE AGN mon- 
itoring data (e.g., Edelson & Nandra 1999). Unmodelled 
variations in the instrument background arc usually less 
than 2 percent of the total observed (sky plus instrument) 
background at energies less than 10 keV (Jahoda et al., 
in prep.''). Ignoring the statistical uncertainty (there was 
adequate signal-to-noise in all observations), a systematic 
uncertainty of <2 per cent should thus be kept in mind for 
all fluxes. Counts were extracted only from the topmost 
PCU layer to maximize the signal-to-noise ratio. All of the 
targets were faint (< 40 ct s~^ PCU~^), so the applicable 
'L7-240' background models were used. Because the PCU 
gain settings changed three times since launch, the count 
rates were rescaled to a common gain epoch (gain epoch 3) 
by calibrating with several public archive Cas A and Crab 
observations. Light curves binned to 16 s were generated 
for all targets over the 2-12 keV bandpass, where the PCA 
is most sensitive and the systematic errors and background 
are best quantified. Light curves were also generated for 
the 2-4 and 7-12 keV sub-bands. The data were binned on 
the orbital time scale; orbits with less than ten 16-second 
bins were rejected. Errors on each point were obtained 
from the standard deviations of the data in each orbital 
bin. Further details of RXTE data reduction can be found 
in e.g., Edelson & Nandra (1999). 

2.2. RXTE sampling 

The observational approach of this project was to quan- 
tify the continuum variability properties of Seyfert 1 galax- 
ies on multiple time scales. This required assembling 
samples that were, to the greatest degree possible, uni- 
formly monitored for proper comparison between sources. 
Sources with a weighted mean count rate significantly be- 
low 1 ct s~^ PCU~^ over the full 2-12 keV bandpass were 
rejected to minimize the risk of contamination from faint 
sources in the field-of-view, to ensure adequate signal-to- 
noise, and to minimize the influence of systematic varia- 
tions in the modeled X-ray background. 

The sampling of the publicly available data was highly 
uneven in general. The original observations were made 
with a wide variety of science goals, leading to a variety of 
sampling patterns and durations. This required us to clip 
light curves to common durations and resample at similar 
rates in order to produce samples with homogeneous sam- 
pling characteristics. For each total light curve, optimum 
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windows of 1 d, 6 d, 36 d, 216 d, and 1296 d (evenly-spaced 
in the logarithm by a factor of 6) were selected. Given 
the original sampling patterns, these windows represented 
a reasonable spread in temporal frequency coverage, and 
yielded a reasonably-sized sample on each time scale. For 
each time scale, light curves shorter than the optimum 
window were rejected. Light curves with long gaps (>l/3 
of the total duration) within the window were also rejected. 
Such gaps reduce the statistical significance of parameters 
derived over the full duration, and interpolating a(;ross 
such large gaps would result in an underestimate of the 
true variability amplitude. For each source, as many us- 
able light curves as possible on each of the five time scales 
were selected from the total light curves. In NGC 3227, 
there was a significant hardening of the spectrum during 
approximately MJD 51900-52000, consistent with a tem- 
porary increase in cold absorption due to a dense cloud 
passing along the line of sight (Lamer, Uttley & McHardy 
2003b); these data were excluded. 

To extract light curves that were sampled as uniformly 
as possible, the light curves were resampled on each of the 
five time scales with a common, optimized rate. This was 
done using an algorithm which selected only data points 
in the original light curve that were separated as close as 
possible to the resampling rate ATgamp, where ATgamp was 
5760 sec (1 satellite orbit), 0.27 d (4 satellite orbits), 1.6 
d, 5.3 d and 34.4 d for the 1, 6, 36, 216, and 1296 d light 
curves, respectively. Starting with the first data point ob- 
served (at time ti), the algorithm selected the data point 
observed closest to times h= ti + ATgamp, ••• , tjv = 
tN-i + ATgamp, where N is the number of points in the 
final, resampled light curve. For light curves that were 
observed with overlapping sampling patterns, portions of 
intense monitoring were not treated differently from the 
rest of the light curve. That is, the algorithm did not 
do any averaging (in bins of size ATgamp) during times 
of intensive monitoring, as that would yield reduced vari- 
ability during that period only relative to the rest of the 
light curve. This allowed the entire light curve to sample 
the variability in a uniform fashion. Resampling at rates 
longer than ATgamp would have resulted in too few points 
in each final light curve, while resampling at significantly 
more frequent rates would have resulted in light curves 
that were not sufficiently uniform, given the original range 
of observing patterns. The final light curves were also re- 
quired to contain at least ^20 points (^15 on the 1 d time 
scale) in order to obtain an accurate estimate of the vari- 
ability amplitude as quantified below; those light curves 
with fewer points were discarded. Light curves with poor 
signal to noise (i.e., due to mean count rates significantly 
less than 1.0) were discarded. Given that many sources 
were observed with overlapping sampling patterns, the fi- 
nal light curves for a given source often share data points 
on multiple time scales and are not completely indepen- 
dent. 

This reduction yielded a total of 27 sources with sam- 
pling on each at least one of the five RXTE time scales. 
This included 86 observations of 18 sources on the 1 d time 
scale, 68 observations of 12 sources on the 6 d time scale, 
19 observations of 12 sources on the 36 d time scale, 78 
observations of 19 sources on the 216 d time scale, and 12 
observations of 9 sources on the 1296 d time scale. Figure 



1 shows the full 2-12 keV RXTE light curves for all 27 
sources, before resampling, and showing the boundaries 
of the sampling windows. Table 1 lists source observa- 
tion parameters, including 2-12 keV luminosity T2-12 and 
black hole mass estimate Mbh, and sampling parameters. 
All source luminosities were calculated using the global 
mean RXTE count rate and using the HEASARC's on- 
line WebPIMMS v.3.4 flux converter assuming an intrin- 
sic power-law with a photon index obtained from either 
previously published spectral fits (e.g., Kaspi et al. 2001, 
Pounds et al. 2003) or the online TARTARUS database 
of ASCA AGN observations (e.g., Nandra et al. 1997a; 
Turner et al. 1999). Luminosities were calculated assuming 
Ho = 70 km s~^ Mpc~^ and Qo = 0.5. All black hole mass 
estimates are reverberation-mapped masses from Kaspi 
et al. (2000) and Wandel, Peterson & Malkan (1999) ex- 
cept NGC 4051, from Shemmer et al. (2003), NGC 3783, 
from Onken & Peterson (2002), NGC 4593, NGC 3516 & 
NGC 3227, from Onken et al. (2003), and Mkn 279, from 
Wandel (2002) and Santos Lleo et al. (2001). Mass esti- 
mates for Ark 564, Mkn 766, MCG-6-30-15, MCG-2-58-22 
are from Bian & Zhao (2003) and the mass estimate for 
PKS 0558-504 is from Wang et al. (2001); these latter two 
works use the empirical Kaspi et al. (2000) relation be- 
tween optical luminosity and BLR size. No reliable mass 
estimate exists for 3C 111 or IRAS 18325-5926. 

2.3. ASCA data 

Short-term ASCA 2-10 keV light curves were obtained 
from the TARTARUS database for the sources with RXTE 
data. The count rates in the light curves provided had 
been combined and averaged between ASCA's two Solid- 
state Imaging Spectrometers (SIS; Burke et al. 1994, Gen- 
dreau 1995) and binned to 16 s. For each source, all avail- 
able light curves longer than 1 d in duration were selected 
from the database; otherwise the longest light curve >60 
ksec in duration was used. The light curves were binned 
on orbital time scales, yielding 51 light curves of 11-15 
consecutive orbital bins for 21 sources. Background light 
curves were similarly binned and subtracted to produce 
net count rate light curves. Table 2 lists source observa- 
tion and sampling parameters for the ASCA data. 



3. ANALYSIS 



3.1. 



Quantifying variability amplitudes 

Fractional variability amplitudes {Fyar] e.g., Vaughan 
et al. 2003b, Edelson et al. 2002) were measured for each 
light curve to quantify the intrinsic variability amplitude 
relative to the mean count rate and in excess of the mea- 
surement noise; 
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un.r — 



(X)2 



(1) 



where S"^ is the total variance of the light curve, {a'^^r) 
the mean error squared and {X) is the mean count rate of 
N total points. The error on F^ar is 
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as discussed in Vaughan et al. (2003b); this error formu- 
lation estimates cFp based on random errors in the data 
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itself, and not due to random variations associated with 
red- noise processes. 

In any red-noise stochastic process there will be random 
scatter in independent estimates of the variance or Fyar 
over multiple realizations of the process. This is a form 
of "weakly non-stationary" behavior inherent in red- noise 
variability processes. Herein, we adopt the definition of 
weak non-stationarity as a description of a variability pro- 
cess whose mean and variance show scatter over multiple 
realizations, but whose underlying PSD remains constant 
over time, with expectation values of Fyar, {Fvar)i remain- 
ing constant over time as well (e.g., Vaughan et al. 2003b). 
In other words, while the expectation value of the square 
of Fyar is equal to the integrated PSD of the underly- 
ing variability process, multiple independent realizations 
of that process will yield a range in estimates of Fyar even 
if the PSD does not change amplitude or shape and {Fyar) 
is constant. Factors of 3 or more in the range of Fyar 
are not uncommon (e.g., Vaughan et al. 2003b). Scatter 
in Fyar is therefore not necessarily indicative of strongly 
non-stationary behavior. For multiple RXTE or ASCA 
light curves for a given source and time scale, the values 
of Fyar were averaged, with the uncertainty on the aver- 
age Fyar determined statistically. However, one needs at 
least 10-20 independent estimates of Fyar to adequately 
test if those estimates are consistent with their expecta- 
tion value (Fyar) (see Vaughan et al. 2003b for detailed 
descriptions of such tests). There arc only three objects 
with enough data for this relatively strong test, NGC 7469, 
IRAS 18325-5926 and MCG-6-30-15 on the 1 d time scales 
(with the RXTE 2-12 keV and ASCA 2-10 keV values con- 
sidered together); in all three cases at least 70% of the in- 
dividual values of Fyar are consistent with {Fyar)- For the 
rest of the sample, when multiple estimates of Fyar were 
made, the measured values were usually reasonably close 
to {Fyar)- This is consistent with weakly non-stationary 
behavior. Thus, these values of Fyar are used hereafter. 

A linear relation between absolute rms variability am- 
plitude and flux has been observed in XRBs (Uttley 
& McHardy 2001) and Seyfert Is (Edelson et al. 2002, 
Vaughan, Fabian & Nandra 2003a). This is a form of 
non-stationary behavior (independent of the weak non- 
stationarity discussed above), as the expectation value of 
the variance is not constant over time. Quantifying vari- 
ability using Fyar removes this trend; since the rms-flux 
relation slope is generally seen to be close to 1, Fyar thus 
would be independent of flux level in the absence of addi- 
tional sources of non-stationarity. 

Table 3 lists Fyar for each RXTE light curve over the 
2-12 keV, 2-4 keV and 7-12 keV bands. Table 4 hsts Fyar 
measured over the 2-10 keV band for the ASCA data. 

3.2. Construction of correlation diagrams 

Figures 2a and 2b displays the values of the logarithm of 
Fyar plotted against ^2-12 and Mbh- Fyar itself will not 
follow a Gaussian distribution, but the logarithm of Fyar 
does follow a distribution that crudely resembles a Gaus- 
sian for red- noise PSD slopes of = sim-1 to -2 (see, e.g.. 
Fig. 8 of Vaughan et al. 2003b). The ASCA data are in- 
cluded and agree well with the 1-d RXTE data; one should 
not expect any significant difference between parameters 
derived over the 2-10 and 2-12 keV bands. Best-fitting 



power-law slopes for each data set were determined using 
the Akritas & Bershady (1996) regression, which accounts 
for measurement errors and intrinsic scatter; the slopes are 
listed in Table 5. The Spearman rank correlation coeffi- 
cients r and probability Pr of obtaining those values of r 
by chance are also listed in Table 5. Also listed in Table 
5 is reduced chi-squared Xr' calculated using the best-fit 
power-law. Because the uncertainties on Fyar were calcu- 
lated in linear space, when calculating x^, the uncertain- 
tics in log space were replaced by the average error for all 
five time scales, 0.047 in the log and 0.044 in the log for 
the Fyar-L2-i2 and Fyar-Msu relations, respectively, to 
avoid unnecessarily weighting towards the longer time 
scales. 

For both sets of relations, the values of r are all nega- 
tive, with the absolute values of r decreasing slightly to- 
wards longer time scales. Such anticorrelations have been 
observed previously in AGNs for ^1 d time scales (Green, 
McHardy & Lehto 1993, Nandra et al. 1997a; O'Neill et al. 
2004). However, the slopes and normalizations of the best- 
fitting logarithmic power law for each data set differ: the 
slopes generally flatten towards longer time scales. The 1- 
d and 216-d time scale Fyar ' L2-12 relations are generally 
consistent with the 1 d and 300 d relations of MEOl. For 
both sets of relations, for all objects, the values of Fyar 
generally increase towards longer time scales, leveling off 
somewhat beyond approximately the 36 d time scale rela- 
tion, but the highest mass and highest luminosity sources 
show the largest increase. Formally, the flts to all of the 
best-fit lines are quite poor, but using the Xr values as 
a measure of intrinsic scatter, the 1 d RXTE Fyar Mbk 
relation shows lower scatter than the corresponding Fyar— 
L2-12 relation, but the scatter is greater in the Fyar Mbh 
relation in each of the remaining five data sets. The sum 
of all six Xr values is also greater for the Fyar-MBu re- 
lations. Within each plot, reduced chi-squared tends to 
decrease towards longer time scales, implying greater in- 
trinsic scatter on the shortest time scales. 

It can be seen from the values of Fyar listed in Table 3 
that most observations (56/68) show stronger variability in 
the 2 4 keV band compared to the 7-12 keV band. For- 
mally, the null hypothesis of the 2-4 keV and 7-12 keV 
excess variances (square of the Fyar) being consistent is 
rejected using an F-test at >90% significance in 12 obser- 
vations and >95% significance in 8 observations. The cor- 
relation diagram for 2-4 keV Fyar {Fyar, so ji) versus 7 12 
keV Fyar {Fyar,hard) for all five RXTE data sets is shown 
in Figure 3. Fyar, so ft and Fyar, hard are well-correlated on 
all time scales. However, it can be seen that the vast ma- 
jority of points lie to the right of the dashed line which 
represents equal variability in the two bands. This shows 
again that most sources exhibit stronger variability in the 
relatively softer band. Values of the Spearman rank cor- 
relation coefficient and probabilities are listed in Table 5. 
Best-fitting power-law slopes were determined using the 
Akritas & Bershady (1996) regression. There is no obvious 
indication that the degree of spectral variability exhibited 
is dependent on the time scale probed. There is no obvious 
scatter trend with time scale, judging from the Xr values. 

Figure 4 shows the ratio of _i£n££/i_ 2-4 keV Fyar I 7- 

var,hard ' 

12 keV Fyar-, plotted against L2-12 for all five RXTE data 
sets. The best-fitting slopes, obtained using the Akritas 
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& Bershady (1996) method, and values of Spearman rank 
r and Pr are listed in Table 5. Also listed in Table 5 are 
the slopes and correlation coefficients for versus 

-T liar, hard 

Mbh (not plotted). The best-fitting slopes are all roughly 
similar; the slopes might be taken as tentative evidence for 
relatively less luminous or less massive sources to be more 
strongly variable in the soft band. There is no obvious 
indication that the degree of spectral variability exhibited 
is dependent on the time scale probed. Judging from the 
Xr values, there is no obvious scatter trend with time scale 
with either or Mbh- The respective sums of the five 

Xr values are approximately equal, implying roughly equal 
scatter in the ^-^2-12 and -Mbh relations. 

Color- flux diagrams, in the which the logarithm of the 
7-12 keV / 2-4 keV count rate hardness ratio (HR) is 
plotted against the logarithm of the geometric mean of 
the count rates in these two bands, are shown in Figure 
5. To minimize the effects of changes in spectral response 
duo to PCA gain epoch changes, only the largest number 
of points within a single gain epoch was selected for each 
source. Light curves of ~300 days in duration, with no re- 
sampling, were used; light curves with less than 70 points 
were discarded. This yielded a sample of 27 light curves 
for 14 sources; date ranges are listed in Table 6. For each 
source, the data were sorted by increasing geometric mean 
and grouped into bins of 16 points; the highest flux bin 
was ignored if it contained less than 10 points. For most 
sources, the data form a continuous, well-defined region. It 
is clear from these diagrams as well that nearly all sources 
soften as they brighten. The two exceptions, which show 
either a slight hardening or no spectral variability with 
flux, are the NLSyl Ark 564 and the radio-quiet quasar 
PG 0804-(-761, as has been reported previously (Edelson 
et al. 2002; Papadakis, Reig & Nandra 2003). Also shown 
in Figure 5 is the best fltting linear fit to the binned data. 
Table 6 lists the mean hardness ratio values (HR) for each 
source. For the sample as a whole, the average of the 27 
mean hardness ratios is 1.06. Ten sources' (HR) values are 
within 20% of the sample average. However, two sources 
are notably softer, the soft-spectrum source Ark 564 and 
the quasar PG 0804-1-761, both of which are usually mea- 
sured to have relatively steep photon indices (e.g., Leighly 
1999; Papadakis, Reig & Nandra 2003). Three sources, 
NGC 3227, NGC 3516, and NGC 4151, are notably harder; 
these soiirces are frequently measured to have relatively 
flat photon indices (e.g., Nandra et al. 1997b; George et 
al. 1998). 

Also listed in Table 6 is a parameter derived directly 
from the slope m of the best linear flt, (j) = 2.0"™, which 
quantifies the decrease in HR for every doubling in geomet- 
ric mean count rate. Multiply-measured values of (HR) 
and (j) for a given object tend to be consistent with each 
other, suggesting that sources do not undergo any radical 
changes in spectral variability behavior over times scales 
of one or two years. is greater than 1 for all sources 
except PG 0804-^761 and Ark 564. It is noted that these 
two sources have the lowest 7-12 keV mean count rates in 
the sample (^^0.3 c/s/PCU). As discussed in more detail in 
the Appendix, it is conceivable that systematic variations 
in the modeled background may contribute greatly to the 
observed variability at such low flux levels, particularly in 



the 7-12 keV band. The RXTE data for PG 0804H-761 
and Ark 564 will therefore not be considered further here. 

Figure 6 shows (j) plotted against L2-12 and Mbh- The 
Spearman rank correlation coefficients are given in Table 
5. The best-fit lines in Figure 6 were again cak;ulated us- 
ing the method of Akritas & Bershady (1996). The best- 
fitting slopes, listed in Table 5, are negative, indicating 
that relatively less luminous or less massive sources display 
a stronger degree of spectral variability per given increase 
in overall fiux. To estimate the scatter, was c;valuatcd 
for each plot and is listed in Table 5; formally, the fits to 
the lines are quite poor, but reduced is slightly lower 
in the (/)-Mbh plot. 

It is noted that other studies (Edelson et al. 2002; Pa- 
padakis, Reig & Nandra 2003) have found Ark 564 and 
PG 0804-1-761 to show hardness ratios that are indepen- 
dent of flux, which would imply values of (f) close to 1. It is 
noted that in both sources, this value would lie reasonably 
close to the observed 0-i2-i2 anticorrelation. Addition- 
ally, the high-mass PG 0804-1-761 would lie close to the 
^Mbh anticorrelation; however. Ark 564 would be a sig- 
nificant outlier if added to the (f)-MBH anticorrelation. 

4. DISCUSSION 

When one uses the fractional variability, Fyar, as a de- 
scription of the intrinsic, underlying variability process, 
certain caveats must be kept in mind when red-noise pro- 
cesses are relevant. Each light curve is an independent re- 
alization of the underlying stochastic process and there will 
be random fluctuations in the measured variance. How- 
ever, in the absence of evidence for strongly non-stationary 
behavior in Seyfert light curves (e.g., §3.1, Markowitz et 
al. 2003a, Vaughan et al. 2003b), it is assumed hereafter 
that the values of Fyar are reasonable quantifications of 
the intrinsic variability amplitude. The reader must keep 
in mind the previously discussed limitations when consid- 
ering such small numbers of F^ar estimates. 

The anticorrelations between variability amplitude and 
source luminosity and between variability amplitude and 
black hole mass seen in previous surveys are confirmed here 
on short time scales. In both sets of anticorrelations, the 
best-fitting power-law slopes gradually decrease towards 
longer time scales. The Fyar values tend to increase to- 
wards longer time scales, however, they tend to saturate 
beyond the 36 d time scale. Consequently, the increase 
in Fyar is greatest for the higher luminosity sources. As 
will be discussed in §4.1, this trend is consistent with a 
scaling of PSD break frequency with some fundamental 
parameter, most likely Mbh- All of the sources exhibit 
stronger variability towards relatively softer energies. Ad- 
ditionally, sub-band Fyar values and color-fiux diagrams 
indicate that less luminous sources have a tendency to 
exhibit more spectral variability overall. These spectral 
variability characteristics are discussed in the context of 
simple X-ray reprocessing models in §4.2. 

4.1. The variability-luminosity-MBH relationship 

Recent PSD studies have yielded PSD breaks on time 
scales of a few days or less; generally, the power-law slopes 
flatten from ^-2 above the break to ~-l below the break. 
Markowitz et al. (2003a) developed a picture in which 
all Seyfert 1 PSDs have the same shape but whose high- 



6 



Seyfert 1 X-ray Variability 



frequency break time scales Tb scale linearly in temporal 
frequency with Mbh- This is consistent with observed 
anticorrelations between Fyar and -Mbh (Papadakis 2004, 
O'Neill et al. 2004) and F^ar and X-ray luminosity (e.g., 
Nikolajuk et al. 2004). Interestingly, though, the PSD 
break frequencies appeared to be less correlated with 2- 
10 keV X-ray luminosity. There are not adequate data to 
construct high dynamic range PSDs for all targets in the 
current sample. However, it is reasonable to assume that 
all Scyfcrts have similar PSD shapes with breaks. Given 
the ranges of luminosity and black hole masses spanned 
by the sample, it is reasonable to assume that the longest 
time scales probes in this survey are exploring variabil- 
ity on temporal frequencies well below the breaks in most 
or all of the sources. This would then explain why the 
variability amplitudes observed tend to saturate at similar 
levels on the longest time scales probed, strongly reducing 
the dependence of F^ar on Mbh or ^2-12. However, the 
data are not able to highly constrain if objects' PSDs con- 
tain a second, low-frequency break, due to the saturation 
of F 

^ var • 

Using the values of Fyar measured here, it is possible to 
further test this picture. By scaling the break frequency of 
a broken power-law model PSD with Mbh, the resulting 
predicted values of F^ar can be compared to the observed 
Fyar values on each of the five time scales to quantitatively 
constrain the Tb A/bh relation. Additionally, scaling the 
break frequency of a model PSD with luminosity can sim- 
ilarly constrain the relation between Tb and L2-12 (the 
bolometric and X-ray luminosities can be related as ap- 
proximately iboi=27-L2-i2; Padovani & Rafanelli 1988). 

For both of these tests, it was assumed that all Seyferts 
have the same singly-broken PSD shape P{f) described 
by P{f) = A{f/f,)-' (for / < /(,), or P(/) = A{f/ft,)-^ 
(for f > fb)- A is the PSD normalization at the high- 
frequency break fb, calculated as 0.01 (Hz~^)//6, a re- 
lation estimated from the A- Mbh and Tb^MnH plots of 
Markowitz et al. (2003a; their figures 12 and 13). A lin- 
ear scaling between Tb and Mbh (or -^2-12) was assumed. 
The Fyar values were calculated by integrating the PSD 
between the temporal frequencies of 1/D (where _D is 1, 
6, 36, 216, or 1296 days) and l/2ATsamp- The values of 
Fyar measured from the observed light curves contain ad- 
ditional contributions to the total variability due to alias- 
ing, which arises from the non-continuous sampling, and 
red-noise leak, which arises due to the presence of vari- 
ability on time scales longer than those sampled. The 
reader is referred to Uttley, McHardy & Papadakis (2002) 
or Markowitz et al. (2003a) for detailed descriptions of 
these distortion effects inherent in PSD measurement and 
variability analysis. The contribution to the total vari- 
ance from aliasing was estimated analytically by integrat- 
ing the model PSD from a frequency of l/(2ATsamp) to 
a frequency of 1/ (2000 s) ; no contribution to the aliased 
power is expected from variations on time scales shorter 
than ~2000 s. Monte Carlo simulations were carried out 
to estimate the contribution to the total variance from red- 
noise leak. For each model PSD, a light curve of length 
SOD, where D is the observed light curve duration, was 
simulated and split into 50 light curves each of length D 
to ensure that variability power from red-noise leak was 
present on the same time scales probed by the observations 



(time scales shorter than D) . The average variance of these 
50 light curves was calculated and compared to the esti- 
mated intrinsic (i.e., no red- noise leak present) variability 
estimated above to estimate the variability contribution 
from red-noise leak. 

In order to more directly study the link between the T^,- 
Mbh and Tb-L2-i2 relations and Fyar, it was necessary to 
remove the influence of the PSD amplitude A at the break 
on Fyar- The accumulation of Seyfert PSDs supports a 
range in the observed values of A (e.g., Uttley et al., in 
prep.). To remove the dependence of the T^^-Mbu and Tb- 
L2-12 relations on A, the ratios Rpvar of values of 2-12 
keV Fyar on six combinations of time scales (i?i/6, denot- 
ing Fyari^ d) / Fyar{Q d), i?i/36, i?l/216, -Rl/1296i -^6/2165 

^6/1296) were considered. The remaining four model ratios 

are all relatively flat across the ranges of Mbh and L2-12 
considered. They do not provide constraints on scaling 
and fitting the model Fyar ratio lines, and are therefore 
excluded from analysis. 

The ratios of predicted Fyar values are plotted as solid 
lines as a function of Mbh and L2-12 in Figures 7a and 
7b, respectively. No arbitrary scaling in the y-direction of 
the resulting values of Rpvar was done. Also plotted are 
the ratios of observed Fyar values; observed values of Fyar 
on the 1 d time scale were combined between the ASCA 
and RXTE data sets by averaging multiple values for each 
source. The predicted Rpvar functions were simultane- 
ously best-fit in the x-direction. The fits indicate that the 
best-fit linear PSD scaling for Figure 7a requires the rela- 
tion Tb (days) = Mbh/W^ '^Mq. The fit is formally quite 
poor, with Xr equal to 57.1 for 48 degrees of freedom. For 
Figure 7b, the linear PSD scaling required is Tb (days) 
= L2-12/ (10^3-^ erg s-i), with Xr equal to 490.2 for 48 
degrees of freedom. The modeling is better overall for 
the Tb-Mnn relation compared to the Tb-T2-i2 relation, 
given the respective values of x^. These two best-fitting 
relations together suggest that the average accretion rate 
for the entire sample is ~5% of the Eddington limit. 

McHardy et al. (2004) suggested that the normaliza- 
tion of a linear Tb A^bh relation may be dependent on 
some other parameter, possibly the accretion rate. Under 
the assumption that the reverberation masses are accu- 
rate, the picture emerging from PSD measurement seems 
to be revealing a bifurcation in Seyfert PSDs. It ap- 
pears that the PSD breaks of some Seyferts lie close to 
a Tb-MBH scaling that is approximately quantified as Tb 
(days) = Mbh/IO^-^Mq (e.g., NGC 3516, NGC 4151, and 
NGC 3783; Markowitz et al. 2003a). This relation extrapo- 
lates 6 7 orders of magnitude to the PSD break of Cyg X-1 
in the low/hard state. Other sources (NGC 4051 and pos- 
sibly other Narrow-Line Seyfert Is, McHardy et al. 2004) 
seem to require a T^-Mbh scaling that is approximately 
Tb (days) = Mbh/IO'^'^ M©. This relation extrapolates 
to the PSD break of Cyg X-1 in the high/soft state, ar- 
guing some connections between these Seyfert s XRBs in 
the high/soft state. The best-fitting linear Tb-Mnn rela- 
tion derived from the present sample lies in between these 
two scalings, though much closer to the low/hard state 
scafing, Tb (days) = Mbb/W^'^ Mq relation. This is con- 
sistent with the idea that the present sample contains a 
mixture of sources from the two groups, but the number 
of sources that scale with Cygnus X-l's high/soft state is a 
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Zdziarski & Grandi 2001). Changes in Tint can also arise 
from changes in optical depth (e.g., Haardt, Maraschi & 
Ghisellini 1997), geometry (e.g., Merloni & Fabian 2001) 
and energy balance (e.g., Zdziarski et al. 2003). However, 
spectral variability studies by McHardy, Papadakis & Utt- 
ley (1998), Shih et al. (2002) and Lamer et al. (2003a) have 
shown that the spectral fit photon index saturates at high 
flux. To explain this effect, McHardy, Papadakis & Utt- 
ley (1998) and Shih et al. (2002) independently proposed 
the "two-component" model c:onsisting of a constant hard 
reflection component superimposed upon a soft coronal 
component that is variable in normalization but constant 
in spectral shape. That is, Tint is constant due to both the 
disk seed and coronal fluxes increasing. As an example, a 
weak dependence of the variability on energy, as has been 
observed in some Narrow-Line Seyfert I's (e.g., Edelson et 
al. 2002), is possible in the context of the two-component 
model if the hard component is absent or extremely weak. 

The color-flux diagrams not only show that Seyfert Is 
generally soften as they brighten, they also tentatively 
suggest that there is more spectral variability for a given 
increase in flux for the relatively less luminous, less mas- 
sive, and more variable overall sources. Additional support 
comes from the marginal anticorrelations between the ra- 
tios of the 2-4 keV and 7-12 keV Fyar and luminosity 
(Figure 4) and Mbh- This trend could be due to some 
variable soft component present in the 2-4 keV band but 
not evident at higher energies; this component could be 
more prominent or more variable in the relatively lower 
luminosity objects. Alternatively, the physical parameters 
which ultimately constrain the amount of observed spec- 
tral variability may themselves be more variable in the 
relatively lower luminosity objects. 

Another possible contribution to this effect may arise 
from the energy dependence of the high-frequency PSD 
(e.g., Papadakis & Nandra 2001, Vaughan, Fabian & Nan- 
dra 2003a, McHardy et al. 2004) . At temporal frequencies 
above the break, PSD slopes tend to increase in slope as 
photon energy increases, typically by ~0. 1-0.2 for a dou- 
bling in photon energy. One would then observe a reduc- 
tion in the ratio of soft to hard X-ray variability in more 
massive or luminous sources, since their PSD breaks ap- 
pear at relatively lower temporal frequencies. However, 
simulations show that such an effect is minor. Simulations 
of 300-day light curves using PSD shapes with energy- 
dependent high-frequency slopes (change in slope by 0.2 
between the two bands), energy-dependent normalization 
A (roughly 50% higher in the soft band; e.g., McHardy 
et al. 2004), and a Tt-MsH relation as per above yield a 
reduction in the ratio of soft to hard F-^ar by ~6% over 
the Mbh range of interest. This corresponds to a change 
in (j) of only ~7%, much smaller than the range observed. 

5. CONCLUSIONS 

This paper extends the results of the first long-term X- 
ray variability survey of MEOl to additional sources and 
time scales, including sampling variability on time scales 

® We note that this last analysis step refers to high/soft state PSD sources and not Narrow-Line Seyfert Is because there may not be a 
one-to-one correspondence between classification as a Broad- or Narrow-Line Seyfert and PSD scaling category. For example, given their mass 
estimates, the PSDs of the Broad-Line Seyfert Is NGC 3227 and Ark 120 are more consistent with scaling with the high/soft state of Cygnus 
X-1 (Uttley et al., in prep.; Marshall et al. 2003), while the PSD of the Narrow-Line Seyfert 1 Ark 564 is more consistent with scaling with the 
low/hard state of Cygnus X-1. This is why NGC 3227 and Ark 120 were included in the high/soft PSD scaling category above and Ark 564 
was not. 



small fraction of the whole sample. Ignoring the five known 
high/soft state PSD sources^ (NGC 4051; MCG 6-30-15; 
NGC 3227, Uttley et al., in prep.; Ark 120, MarshaU et 
al. 2003; PG 08044-761, Papadakis, Reig & Nandra 2003) 
indeed gives a slightly lower scaling constant in the best- 
fltting Tb-MBH relation, (days) = Mbk/W'^-'^Mq {Xr 
= 59.4 for 28 degrees of freedom). 

As mentioned previously, the large amount of scatter in- 
herent in Fvar complicates the present analysis. Estimates 
of Fya,- for a given source will contain scatter even when 
{Fyar) is constant, due to the stochastic nature of red- noise 
variability processes. Moreover, not all Seyfert PSDs are 
exactly identical in PSD shape and amplitude, meaning 
that there will be some scatter in (Fyar) from one object 
to the next, even when Fyar is consistently measured over 
identical sampling windows. For instance, fixing the high- 
frequency PSD slope and break frequency while doubling 
the PSD amplitude A will increase (Fyar) by 41 percent for 
all time scales studied. For fixed A and break frequency 
fixed at 10^^ Hz, steepening the high-frequency PSD slope 
from -2.0 to -2.5 will decrease (Fyar) by 1.4 and 2.2 (de- 
creases of 0.16 and 0.34 in the log) on the 6 d and 1 d time 
scales, respectively. Finally, the aforementioned bifurca- 
tion in PSD break frequencies, corresponding to scaling 
with either the high/soft state or low/hard state of Cyg 
X-1, introduces scatter in {Fyar) on time scales longer than 
Tb. For time scales of a year or more, for the black hole 
masses and PSD break frequencies of interest, {Fyar) will 
change by a factor of ~3 or more. There hence is intrinsic 
scatter in {Fyar) at both long and short time scales due 
to these effects. Assuming that the values of Xr given in 
Table 5 are an adequate characterization of the intrinsic 
scatter in the Fyar-L2-i2 and Fyar^Msn relations, one 
could speculate that the increased scatter towards shorter 
time scales may indicate that the range of high-frequency 
PSD slopes contributes more to the overall scatter than 
the low-frequency PSD bifurcation. However, removal of 
the five known high/soft state PSD sources fails to reduce 
scatter at long time scales, and it remains difficult to iden- 
tify the dominant source of intrinsic scatter. 

4.2. Spectral variability 

The majority of the Seyferts sampled show stronger vari- 
ability towards softer energies, as seen from a comparison 
of the 2-4 keV and 7-12 keV Fyar values, and from the 
color-flux diagrams. Such behavior is consistent with the 
well-documented property of Seyfert Is to soften as they 
brighten. Some works have suggested spectral pivoting of 
the coronal power law about some energy above 10 keV 
as the explanation for Seyferts' softening as they brighten 
(e.g., Papadakis et al. 2002). Thermal Comptonization 
models predict changes in the intrinsic spectral slope of 
the coronal component. Tint- In the case of a coronal cloud 
that is fed by a variable soft photon seed flux, held at con- 
stant optical depth, and not pair-dominated, an increase in 
seed flux will lower the electron temperature of the corona 
and steepen the X-ray spectrum (e.g., Maraschi et al. 1991, 
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well below the putative PSD breaks in Scyferts. The 
well-studied luminosity-variability amplitude anticorrela- 
tion and the anticorrelation between black hole mass and 
variability amplitude are confirmed on short time scales. 
Variability amplitudes increase towards longer time scales, 
consistent with red-noise variability, but the relatively 
more luminous and more massive sources show the great- 
est increase. For both sets of anticorrelations, the best- 
fitting slopes decrease towards longer time scales. These 
trends are consistent with a simple scaling of PSD break 
frequency with black hole mass as suggested by Markowitz 
et al. (2003a) and McHardy et al. (2004), with F^ar sat- 
urating on time scales below the PSD breaks. The best- 
fitting time scale-mass relation is quantified as Tb (days) = 
Mbh/IO^'^Mq, and the best-fitting time scale-luminosity 
relation is quantified as Tb (days) = L2-i2/(10'^^'^ erg 
s~^), implying an average accretion rate for the entire sam- 
ple of ^5% of the Eddington limit. The measurement of 
a larger number of Seyfert PSDs at low temporal frequen- 
cies and additional accumulation of Fyar measurements on 
multiple time scales for a given object will further clarify 
the relations between PSD break time scale, PSD normal- 



ization, Fyan black hole mass and luminosity. 

Nearly all the observations show relatively stronger vari- 
ability towards softer energies, as seen from the values of 
F„ar- Color-flux diagrams additionally show that sources 
soften as they brighten. The color-flux diagrams also ten- 
tatively suggest that sources with relatively lower lumi- 
nosities or black hole masses display a larger range of spec- 
tral variability for a given increase in total X-ray flux. 
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APPENDIX 

INFLUENCE OF PCA BACKGROUND MODELING ON MEASURED VARIABILITY PROPERTIES 

Because RXTE is a non-imaging instrument, the background must be modeled. However, for very faint targets, 
including most soft-spectrum Seyferts, the estimated PCA background count rate is greater than the source count rate. 
Small systematic errors in the background model will thus cause proportionally larger problems for soft-spectrum and low 
count rate sources. In extreme cases, uncertainty in the background model can lead to incorrect characterization of the 
true variability (e.g., the iiATi? observation of the soft-spectrum source TON S180, Edelson et al. 2002). This Appendix 
explores the influence of the background subtraction by examining the measured variability characteristics as a function 
of source count rate. 

Figure 8 shows the logarithm of 2-12 keV F^ar plotted against the logarithm of the count rate for all sources and 
time scales (2-10 keV for the ASCA data). For multiply-observed sources on each time scale, the values of F^ar and 
count rate obtained before averaging were used in order to explore the widest range of count rates possible. Spearman 
rank correlation coefficients and probabilities are listed in Table 7. Weak to moderate anticorrelations are evident on all 
time scales. These are not the result of any correlation between coimt rate and luminosity. Spearman rank correlation 
coefficients and probabilities are listed in Table 7 for source mean count rate as a function of both L2-12 and A'/bh! m 
general, count rate is seen to be uncorrelated with either source parameter. The anticorrelations are, however, the result 
of the inclusion of the narrow-line/ soft-spectrum sources. RXTE, lacking coverage below 2 keV, generally cannot observe 
most soft-spectrum sources. Invariably, this class of objects will yield lower 2-12 keV count rates compared to normal, 
broad-line Seyferts. However, these objects also tend to be more variable than broad- line Seyferts (e.g.. Turner et al. 
1999). Recalculation of the Spearman rank correlations, excluding the six narrow- line/soft spectrum sources, shows the 
above anticorrelations to be substantially weakened on most time scales. However, on the 6 and 36 d time scales, it is 
necessary to additionally exclude NGC 3227. The 2000 intensive monitoring campaign of NGC 3227, from which the 6 d 
and 36 d light curves are derived, happened to catch this highly variable source in a relatively low flux state (as shown 
in Figure 1). Removal of these data points further weakens the anticorrelations in those two plots. Overall, there is no 
evidence that the measured variability characteristics of low count rate sources are affected by the RXTE background 
modeling on any time scale. 

Figure 9 shows the ratio of 2-4 keV F^ar I 7-12 keV Fyar is plotted against the geometric mean of the count rates 
in these two bands for all five RXTE time scales. Again, for multiply-observed sources on each time scale, the values of 
Fvar.softi Fyar.hard and count rate obtained before averaging were used. The Spearman rank correlation coefficients and 
probabilities are listed in Table 7. The ratio is seen to be generally independent of count rate for all time scales. 

The parameter (p is plotted against average geometric mean count rate in Figure 10, and seen to be independent of 
count rate for most sources. However, Ark 564 and PG 0804-1-761, denoted by open circles in the figure, have values of (f> 
less than 1. Systematic background errors may be biasing the estimate of the hardness ratio; this may be an artifact of 
the low count rates for these sources, especially in the hard band (7-12 keV count rates per PCU are about 0.3 counts 
sec~^ for both sources). Spearman rank coefficients and probabilities are given in Table 7 with and without Ark 564 and 
PG 0804-1-761. 
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Seyfert 1 X-ray Variability 



Table 1 



Source and RXTE sampling parameters 



Time 


Source 




log(L2_12) 


log(MBH ) 


MJD 


Num. 


Mean 


Mean 


Scale 


Name 


z 


(erg s ^) 


(M©) 


Range 


Pts. 


— 1 

c s 


S/N 


1 d 


PKS n^'iS-^n4 


0.137 


44.90 


7.7 


50734 31-50735 24 


14 


1.29 


48 




3C 111 


0.049 


44.63 




50530 Q3-50531 66 


w 


4.36 


102 












51982 37-51983 43 


13 


7.95 


72 












51983 43-51984 49 


12 


7.94 


76 




Mkn 509 


0.034 


44.33 


9 2±1 1 


52012.20-52013.26 


15 


6.09 


75 




3C 120 


0.033 


44.24 


7.481^2^ 


52621.91-52622.90 


12 


5.47 


69 










52622.90-52623.89 


13 


5.61 


66 












52623.89-52624.88 


13 


5.98 


79 












52624.88-52625.93 


11 


5.70 


82 












52625.94-52626.99 


11 


5.29 


87 












52628.31-52629.10 


11 


4.81 


68 












52629.37-52630.35 


12 


6.54 


83 












52678.12-52679.15 


11 


6.01 


82 












52679.18-52681.22 


14 


5.73 


73 




MCG-2-58-22 


0.047 


44.24 


8.5 


50797.97-50798.91 


15 


3.59 


100 












51486.10-51487.10 


12 


3.19 


69 




Ark 120 


0.032 


44.02 


8.27l°:?? 


51163.37-51164.37 


12 


3.94 


101 










51164.71-51165.51 


13 


3.83 


106 




Mkn 279 


0.030 


43.82 


7.4 


52412.29-52413.21 


14 


1.40 


18 




NGC 5548 


0.017 


43.63 


J 07+0-08 
— 0.07 


50984.88-50985.88 


14 


9.00 


183 










50985 89-50986 88 


12 


9.15 


202 




Ark 564 


0.025 


43.62 


6.1 


50440.66-50441 .66 


15 


2.21 


53 




IN VjL^ / 4oy 


U.Ulo 


40.00 








0. io 


DO 












5Uz45.U4— 5uz4d.u4 


14 


3.31 


62 












50246.05-50247.04 


15 


2.75 


57 












50247.05-50248.04 


14 


3.38 


60 












50248.05-50249.05 


15 


3.48 


57 












50249.05-50250.05 


14 


3.65 


64 












50250.06-50251.05 


13 


3.41 


54 












50251.06-50252.05 


15 


3.26 


55 












50252.06-50253.06 


13 


2.38 


48 












50253.06-50254.06 


13 


2.73 


53 












50254. 07-50255.06 


13 


3.09 


53 












o02d5.0 (^— o02od.0d 


13 


2.56 


44 














w 


3.42 


58 












50269.44-50270.43 


11 


3.40 


68 












50270.64-50271.58 


11 


3.03 


67 












50271.64-50272.58 


11 


2.40 


53 












50272.65-50273.58 


13 


2.74 


55 












50273.65-50274.57 


13 


2.96 


65 












50274.65-50275.65 


14 


3.08 


68 




IRAS 18325-5926 


0.020 


43.36 





50807.10-50808.09 


15 


1.67 


47 












50808.10-50809.10 


15 


2.48 


67 












50865.71-50866.70 


15 


1.80 


58 












50866.70-50867.70 


14 


1.61 


59 




NGC 3783 


0.010 


43.31 


6 

— 06 


50112.73-50113.59 


13 


5.84 


100 




Mkn 766 


0.013 


43.27 


5.9 ' 


52036.02-52037.08 


12 


4.01 


43 




NGC 4593 


0.009 


43.03 


6.82+2?5 

— 0.57 


52089.03-52089.96 


14 


4.58 


73 










52455.91-52456.96 


12 


5.18 


60 






n Ann 

u.uuy 


4o.UU 


7 r,o+0.07 

' •^'^-0.10 


ouooy.yo- ouoyu.y / 


10 


A /in 
4.4U 












50590.98-50591.97 


15 


4.18 


125 












50591.98-50592.98 


15 


4.40 


145 












50592.98-50593.98 


15 


3.91 


136 












50916.30-50917.30 


15 


5.90 


151 












50917.30-50918.30 


15 


5.74 


141 












50918.31-50919.30 


15 


5.39 


143 












52009.36-52010.35 


15 


3.23 


64 




MCG-6-30-15 


0.008 


42.99 


6.2 


50664.14-50665.14 


15 


4.33 


106 












50665.15-50666.15 


15 


4.47 


104 












50666.15-50667.15 


15 


4.50 


101 












50667.16-50668.15 


15 


5.08 


90 












50668.16-50669.15 


15 


4.38 


94 












50669.16-50670.16 


15 


4.04 


94 












50670.16-50671.16 


15 


4.00 


77 












50671.16-50672.16 


15 


5.01 


99 












51378.13-51379.12 


15 


4.69 


61 












51379.13-51380.12 


15 


5.41 


72 
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Time 


Source 




log{L2-12) 


log(AfBH) 


MJD 


Num. 


Mean 


Mean 


Scale 


Name 


z 


(erg s~^) 


(Mo) 


Range 


Pts. 


C S""*" 


S/N 



NGC 3227 0.004 42.20 



NGC 4051 0.002 41.44 



7.56; 



+0.06 
0.08 



c; 7+0.3 



6 d 3C 120 
Fairall 9 



IC 4329a 



0.033 44.24 
0.047 44.17 



7.48 
7.92 



+0.21 
0.28 

+0.11 
0.32 



Mkn 279 
NGC 5548 



0.016 



0.030 
0.017 



43.99 



43.82 
43.63 



6.85 



+0.65 
6.85 



Ark 564 0.025 43.62 



NGC 7469 0.016 43.38 



NGC 3783 0.010 43.31 



MCG-6-30-15 0.008 42.99 



7.4 

7 Q7+0.08 
'•»'_0.07 



6.1 



6.88 



+0.29 
6.88 



6.94 



+0.05 
0.06 



6.2 



NGC 4151 0.003 42.59 



NGC 3227 0.004 42.20 



7.08; 



+0.23 
0.38 



7.56; 



+0.06 
0.08 



51380.13-51381.12 


15 


5.56 


66 


51381.13-51382.12 


14 


6.74 


77 


51382.13-51383.12 


14 


5.30 


77 


51383.13-51384.12 


15 


5.48 


68 


51384.13-51385.12 


15 


5.86 


70 


51385.13-51386.13 


15 


5.61 


75 


51386. 14-51387. 13 


15 


4.69 


68 


51387.13-51388.13 


15 


5.58 


68 


51736.06-51737.06 


13 


4.37 


42 


51930.73-51931.58 


13 


5.55 


76 


50405.69-50406.68 


12 


6.05 


135 


50406.69-50407.69 


15 


5.68 


129 


50407.69-50408.69 


15 


5.44 


114 


50408.76-50409.56 


12 


5.32 


130 


50430.43-50431.42 


14 


0.67 


17 


50431.82-50432.82 


15 


1.72 


33 


51626.99-51627.99 


14 


1.81 


28 


52044.39-52045.38 


15 


2.11 


27 


52046.51-52047.50 


15 


2.11 


28 


52621.94-52627.95 


24 


5.57 


71 


52144.89-52150.93 


22 


2.54 


20 


52151.00-52157.04 


22 


2.44 


17 


52157.17-52163.21 


21 


2.75 


26 


52163.34-52169.45 


23 


2.22 


22 


52169.58-52175.62 


23 


2.07 


23 


52830.79-52836.65 


27 


13.08 


72 


52837.05-52843.17 


27 


11.46 


62 


52843.37-52849.30 


30 


12.33 


63 


52849.50-52855.75 


26 


12.75 


65 


52855.95-52862.08 


24 


12.79 


65 


52412.32-52418.07 


18 


1.44 


18 


52091.66-52097.70 


22 


5.95 


48 


52097.77-52103.87 


23 


6.26 


51 


52104.01-52109.98 


23 


5.30 


45 


52110.11-52116.22 


23 


4.55 


37 


52116.35-52122.46 


22 


3.19 


26 


51694.85-51700.85 


23 


2.27 


20 


51700.98-51707.04 


22 


2.02 


17 


51707.17-51713.50 


25 


2.20 


18 


51713.50-51719.56 


22 


2.22 


21 


51719.69-51725.68 


21 


2.26 


20 


50244.07-50250.09 


24 


3.26 


59 


50250 15-50256 17 


25 


2.92 


45 


50256 30-50262 38 


25 


2.65 


60 


50262.45-50268.46 


24 


3.41 


88 


50268.53-50274.55 


25 


3.01 


71 


51960.17-51966.28 


24 


6.13 


43 


51966.41-51972.53 


23 


6.27 


46 
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21 


6.44 


43 


50664 1 8-50670 1 9 


24 


4.46 


90 


51378.16-51384.16 


24 


5.45 


70 


51622.71-51628.97 


23 


5.13 


58 


51629.43-51635.43 


18 


4.26 


50 


51635.70-51641.69 


23 


4.74 


57 


51642.22-51648.22 


22 


4.96 


53 


51648.42-51654.41 


20 


4.07 


43 


51654.68-51660.61 


20 


4.65 


39 


51661.27-51667.27 


22 


4.26 


40 


51667.40-51673.40 


18 


4.98 


51 


51673.60-51679.66 


21 


4.91 


44 


51679.86-51685.78 


21 


4.52 


35 


51870.64-51876.69 


24 


8.26 


43 


51876.75-51882.74 


21 


8.87 


52 


51882.80-51888.86 


22 


8.93 


63 


51889.12-51895.10 


22 


7.12 


53 


51895.24-51901.22 


23 


7.92 


64 


51636.59-51642.51 


22 


3.69 


43 
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log(MBH) 
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z 


(erg s~^) 
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51642.85-51648.91 


23 


3.07 


39 








51649.11-51655.04 


21 


2.58 


33 








51655.30-51661.30 


17 


0.99 


14 








51661.56-51667.63 


23 


0.89 


11 








51667.96-51674.09 


20 


1.35 


17 








51674.62-51680.81 


19 


1.34 


14 








51681.01-51687.07 


23 


2.52 


23 








51687.34-51693.33 


22 


3.06 


27 








51693.60 51699.60 


21 


3.49 


33 


0.002 


41.44 


c 7+0.3 
°- ' -0.4 


51665.34 51671.34 


24 


2.25 


27 






51671.54-51677.60 


19 


2.09 


24 








51678.06-51684.06 


24 


2.82 


24 








t^lRHA 19 t^lRQfl 

1^ oioyu.uo 


9*^ 


A 90 


ox 








51690.38-51696.64 


24 


3.76 


32 








51696.91-51702.90 


23 


2.89 


24 








51703.10-51709.36 


26 


3.04 


23 








51709.63-51716.55 


23 


3.23 


24 








51717.15-51723.27 


23 


2.36 


20 








51723.34-51729.33 


18 


2.38 


17 








0ZU4Z.OD— 0ZU4c5.o4 


1 Q 

iy 


9 1 Q 


97 


u.uou 


44 48 




50990 fi'? 50957 05 


9'^ 


2 72 


70 


0.033 


44.24 


T ,10+0.21 


50458.53-50494.49 


22 


5.80 


107 


0.047 


44.17 




52144.89-52179.00 


22 


2.. 38 


20 


0.016 


43.99 




50665.82-50701.87 


22 


13.14 


167 






52830.79-52864.84 


22 


12.61 


63 


0.017 


43.63 


7.97in? 


52091.66-52125.44 


22 


4.97 


42 


0.025 


43.62 


6.1 


51694.85-51726.48 


20 


2.29 


21 


0.016 


43.. 38 


6.88«-f, 


50244.07-50276.02 


21 


3.04 


62 


0.009 


43.00 


7 9.,+0.07 
' •^•^-0.10 


50523.03-50556.60 


21 


6.18 


85 






50557.13-50590.94 


22 


4.43 


65 








50591.01-50624.88 


25 


3.90 


51 








50627.04-50659.11 


21 


4.08 


52 


0.008 


42.99 


6.2 


51622.71-51655.68 


21 


4.50 


47 








51655.95-51688.58 


20 


4.49 


38 


0.003 


42.59 


7.08t;];l 


51870.64-51904.88 


22 


8.84 


58 


0.004 


42.20 


7-56l°:°| 


51636.59-51669.36 


22 


2.16 


27 






51669.62-51702.59 


21 


2.64 


24 


0.002 


41.44 


5 7+°-3 


51665.34-51698.11 


21 


3.05 


29 






51698.24-51731.00 


21 


2.68 


19 


0.1 


44.59 


8.21 ± 0.04 


51610.61-51826.75 


39 


1.28 


17 


0.056 


44.48 


°-'-"-0.48 


51 1 86 05-51 402 29 


38 


2.24 


34 








oo 


'\ QQ 

o.yy 


Oo 










O ( 


QS 

o.yo 


^0 

ou 


0.033 


44.24 


-7 ,|Q + 0.21 


50812.09-51034.89 


23 


4.46 


59 






51039.95-51256.43 


32 


4.89 


65 








51260.23-51479.30 


31 


5.31 


55 








52334.94-52550.15 


32 


4.82 


37 


0.047 


44.17 


7 92+° " 
' -^^-0.32 


50390.63-50598 . 35 


36 


2.81 


56 






50604.01-50808.01 


36 


2.42 


47 








Klicn 5Q— 51 ^Q^ 7s 
oj.iou.uy oj-oyo. / o 




1 7Q 
1. / y 


9fi 








51398.20-51611.51 


49 


1.42 


21 








51615.81-51829.08 


46 


1.33 


15 








51833.52-52046.64 


50 


1.49 


16 








52050.96-52264.31 


44 


2.51 


23 








52268.55-52481.78 


46 


2.56 


26 








52486.19-52699.54 


42 


2.01 


19 


0.032 


44.02 


Q 07+0.09 
8-^'-0.11 


51026.23-51242.33 


39 


3.17 


48 






51425.10-51644.10 


40 


3.28 


40 


0.016 


43.99 


6.85«;i 


52831.32-53048.37 


38 


13.72 


78 


0.035 


43.88 


6.89«-i? 


51610.59-51830.22 


22 


2.02 


22 


0.026 


43.83 




51684.46-51905.06 


26 


3.42 


28 


0.017 


43.63 


7 Q7+0.08 
'-0.07 


50208.07-50426.16 


29 


6.12 


70 






50437.54-50649.39 


17 


4.30 


70 








50661.14-50871.76 


15 


5.60 


83 








50886.70-51095.11 


20 


6.87 


100 








51110.70-51328.96 


40 


6.39 


73 








51333.20-51550.72 


48 


5.99 


63 



36 d 



NGC 4051 



3C 390.3 
3C 120 
Fairall 9 
IC 4329a 

NGC 5548 
Ark 564 

NGC 7469 
NGC 3516 



MCG-6-30-15 0.008 

NGC 4151 
NGC 3227 



NGC 4051 

216 d PC 0804+761 
3C 390.3 



3C 120 



Fairall 9 



Ark 120 

IC 4329a 
Mkn 110 
Mkn 590 
NGC 5548 



Markowitz & Edelson 
Table 1 — Continued 



Time 


Source 




log(L2-12) 


log(MBH) 
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z 
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51554.98- 


-51772.72 


48 


3.66 


40 










51776.90- 


-51994.44 


46 


4.24 


39 










51998.80- 


-52216.32 


45 


3.39 


32 












c;0/i QQ on 


4o 


Q OQ 
O.ZO 


Qn 
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52442.49- 


-52659.94 


47 


1.69 


18 


Ark 564 


0.025 


43.62 


6.1 


51179.58- 


-51392.92 


43 


1.85 


26 










51397.23- 


-51610.51 


44 


1.86 


23 










51614.75- 


-51827.95 


41 


1.76 


16 










51832.24- 


-52045.71 


47 


1.88 


17 










52050.01- 


-52263.32 


45 


1.90 


16 










52267.61- 


-52485.13 


48 


1.74 


16 










52489.41- 


-52702.88 


46 


1.78 


18 


Mkn 335 


0.026 


43.62 


6.58+°" 

— 0.13 


51661.98- 


-51880.63 


23 


2.14 


21 


ivlKn ( y 


n noo 
yj.xjzz 


A'i Afl 
40.4D 


°-^J^-0.35 


olDlU.oy- 




00 

zz 


l.o4 


on 


NGC 7469 


0.016 


43.38 


„ „n+0 29 

6-88i;!:S 


52737.04- 


-52954.51 


43 


2.94 


25 


NGC 3783 


0.010 


43.31 


R 04+0-05 


51180.55- 


-51398.15 


49 


7.89 


82 








51402.34- 


-51615.79 


47 


7.80 


77 










51624.30- 


-51841.91 


45 


8.11 


62 










51846.14- 


-52063.78 


47 


7.73 


56 










52068.06- 


-52289.92 


47 


7.80 


55 


NGC 3516 


0.009 


43.00 


7 90+0.07 
' •■^"J— 0.10 


50523.03- 


-50731.55 


41 


4.52 


57 








50740.07- 


-50949.15 


35 


5.83 


79 










50953.41- 


-51162.41 


46 


4.19 


64 










51166.60- 


-51375.86 


46 


3.74 


50 










51379.80- 


-51593.40 


47 


2.78 


44 


MCG-6-30-15 


0.008 


42.99 


6.2 


50411.95- 


-50622.44 


15 


4.81 


64 










50639.80- 


-50858.59 


17 


5.20 


66 










51870.90- 


-52076.48 


37 


4.89 


38 










52078.47- 


-o22o4.oo 


31 


5.65 


40 










52286.87- 


cro /I r\i OA 

-5z491.o4 


37 


5.43 


42 










52494.07- 


-52701.47 


32 


5.31 


40 


NGC 4151 


0.003 


42.59 




51179.56- 


-51397.20 


50 


18.31 


139 








51401.48- 


-51619.06 


47 


16.49 


121 










51623.29- 


-51840.99 


49 


18.47 


101 


NGC 3227 


0.004 


42.20 


(.00_g Qg 


51180.48- 


-51382.74 


31 


4.92 


60 








51447.16- 


-51663.36 


32 


3.02 


40 










51663.56- 


-51879.93 


27 


2.10 


20 










52174.71- 


-52390.56 


38 


4.19 


35 










52392.54- 


-52609.78 


30 


4.38 


37 


NGC 4051 


0.002 


41.44 


c 7+0.3 
"•'-0.4 


50196.52- 


-50411.78 


26 


2.95 


36 








50427.85- 


-50647.05 


18 


2.02 


31 










50663.15- 


-50888.20 


17 


1.11 


18 










50899.34- 


-51124.11 


17 


1.34 


22 










cr -| -1 Q '7 no 

Olio / .92- 


cr -| '3 r n /I 
-OlODO.Ul 


20 


1.91 


26 










51394.87- 


-51609.53 


15 


1.98 


29 










51611.45- 


-51828.02 


41 


2.17 


19 










51829.15- 


-52045.39 


39 


1.25 


11 










52045.42- 


-52262.34 


38 


2.11 


19 










52264.38- 


-52478.81 


38 


2.51 


23 










52480.99- 


-52698.63 


43 


2.46 


20 


3C 120 


0.033 


44.24 


7.48l°|i 


50458.53- 


-51563.21 


25 


5.07 


71 


Fairall 9 


0.047 


44.17 


' -^^-0.32 


51180.59- 


-52477.68 


39 


1.86 


22 


NGC 5548 


0.017 


43.63 


7 Q7+0.08 
'•y'-0.07 


50208.07- 


-51473.99 


40 


5.66 


71 








51478.23- 


-52749.70 


38 


3.51 


35 


Ark 564 


0.025 


43.62 


6.1 


51179.58- 


-52476.67 


38 


1.86 


18 


NGC 3783 


0.010 


43.31 




51180.55- 


-52375.11 


35 


7.85 


59 


NGC 3516 


0.009 


43.00 


7 90+0.07 
'•^•^-0.10 


50523.03- 


-51819.48 


35 


4.22 


58 


MCG-6-30-15 


0.008 


42.99 


6.2 


50159.80- 


-51423.91 


41 


5.13 


61 










51455.99- 


-52701.47 


37 


5.62 


42 


NGC 3227 


0.004 


42.20 


7 cfj+o.oe 

'•^"-0.08 


51180.48- 


-52476.57 


34 


3.20 


30 








50196.52- 


-51438.88 


39 


1.76 


25 


NGC 4051 


0.002 


41.44 


r; 7+0.3 
'^■'-0.4 


51448.15- 


-52702.59 


38 


1.92 


18 



Note. — The targets arc ranked by 2—12 kcV luminosity (col. [4]). Redshifts (col. [3]) were obtained from the NED database. 
References for black hole mass estimates (usually the reverberation-mapped estimate, e.g., Kaspi et al. 2000 and Wandel, 
Peterson & Malkan 1999) are given in §2.2 of the text. Col. (7) is the number of points in the light curve after clipping and 
resampling to a common rate. Col. (8) is the mean 2-12 keV count rate per PCU. Col. (9) is the signal-to-noise. 
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Table 2 

asca sami'llxc p.,\.r ameters 



Source 


MJD Date 


Sequence 


Num. 


Mean 


Mean 


Name 


Range 


ID Number 


Pts. 


cs-i 


S/N 


3C 111 


50126.54-50127.49 


74087000 


15 


0.74 


34 


PG 0804+761 


50756.55-50757.55 


75058000 


15 


0.22 


21 


3C 390.3 


49307.84 49308.94 


70005000 


15 


0.35 


24 


Mkn 509 


49471.39-49472.39 


71013000 


15 


1.02 


39 


3C 120 


49400.66-49401.66 


71014000 


15 


1.06 


45 


MCG-2-58-22 


49132.89-49133.70 


70004000 


13 


0.24 


12 




50600.47-50601.47 


75049000 


15 


0.71 


33 




50797.98 50798.98 


75049010 


15 


0.70 


33 


Fairall 9 


49688.22-49688.89 


73011000 


11 


0.61 


22 


Ark 120 


49624.79-49625.79 


72000000 


15 


0.65 


31 


IC 4329a 


49214.29-49215.24 


70005000 


15 


1.54 


42 


NGC 5548 


49195.65 49196.60 


70018000 


14 


0.92 


25 




50984.60 50985.60 


76029010 


15 


1.57 


44 




50985.60-50986.60 


76029010 


15 


1.72 


62 




50986.60-50987.60 


76029010 


15 


1.54 


32 


Ark 564 


50440.64-50441.64 


74052000 


15 


0.57 


23 


NGC 7469 


49323.32 49324.06 


71028010 


11 


0.76 


30 


IRAS 18325-5926 


49241.51 49242.51 


70015000 


15 


0.13 


9 




50534.59-50535.59 


75024000 


15 


0.43 


18 




50535.59-50536.59 


75024000 


15 


0.44 


21 




50536.59 50537.59 


75024000 


15 


0.41 


18 




50537.59-50538.59 


75024000 


15 


0.35 


17 




50538.59-50539.59 


75024000 


15 


0.32 


17 


NGC 3783 


50278.28-50278.96 


74054020 


11 


1.55 


41 


Mkn 766 


49339.13-49340.02 


71046000 


14 


0.46 


18 


NGC 4593 


49361.04-49362.04 


71024000 


13 


0.76 


34 


NGC 3516 


49444.15-49445.05 


71007000 


13 


1.57 


25 




50915.94-50916.94 


76028000 


15 


0.87 


36 




50916.94-50917.94 


76028000 


15 


0.91 


42 




50917.94 50918.94 


76028000 


15 


0.89 


37 




50918.94-50919.94 


76028000 


15 


0.72 


38 


MCG-6-30-15 


49177.27-49178.27 


70016000 


14 


0.93 


23 




49199.39-49200.39 


70016010 


15 


0.76 


21 




49556.25-49557.25 


72013000 


15 


1.00 


23 




49557.25-49558.25 


72013000 


15 


1.29 


21 




49558.25-49559.25 


72013000 


15 


1.09 


30 




49559.25-49560.25 


72013000 


15 


0.88 


27 




50663.95-50664.95 


75006000 


14 


0.76 


21 




50664.95 50665.95 


75006000 


15 


0.84 


25 




50665.95-50666.95 


75006000 


15 


0.85 


28 




50667.70-50668.70 


75006010 


15 


0.75 


25 




50668.70-50669.70 


75006010 


15 


0.86 


22 




50669.70-50670.70 


75006010 


14 


0.69 


22 


-\jr~i /I 1 C 1 

INLiO 4iOi 


4yo4 / .14— 4yo4o.UU 


i ouiyuuu 




i.yy 


/ 




49848.00-49848.80 


73019000 


12 


2.42 


69 




49848.80-49849.60 


73019000 


12 


2.52 


66 


NGC 3227 


49115.14-49116.12 


70013000 


15 


0.78 


37 




49852.06-49853.03 


73068000 


13 


0.51 


29 


NGC 4051 


49102.94 49103.89 


70001000 


14 


0.28 


13 




49510.61-49511.54 


72001000 


14 


0.55 


15 




49511.54-49512.47 


72001000 


14 


0.51 


15 



Note. — Targets are ranked by 2-12 keV luminosity, given in Table 1. Col. (4) is the 
number of points in the ASCA light curve after orbitally binning. Col. (5) is the mean 
count rate averaged between both SIS instruments. The SIS data were unusable for the 
observation of NGC 4593 and the first observation of NGC 3227; the GIS data were used 
and GIS count rates were converted to SIS count rates using the online W3PIMMS tool. 
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Table 3 

Derived Variability Parameters for the RXTE light curves 



i ime 


Source 


10 ^^^^\T 

Z—lZ kcv 


2-4 KeV fyar (,/oJ 


(-LZ K.e\ I'yar (/oj 








Scale 


Name 


Fvar (%) 


(-^uar,5o/t ) 


{,-^var,hard) 


^var,hard 


1 H 

1 CI 




191 -1- f! 
iz . 1 zc u.u 


11 c _|_ n Q 

11.0 HZ \j.iJ 


1 Q -U 1 K 
10. y ZIZ 1.0 


0.85 


± 


0.11 




"^P 111 
ov^ J. J. J. 


i.x HI u. / 


4-1-99 


1 .'± ZIZ 1 .'-± 


0.31 


± 


1.58 




1V1K.11 0\}tJ 


1.0 HZ u.^ 


O.U zc u.o 












izu 




4 1 _i_ fl Q 
4:. 1 ziz U.y 


4 c: _|_ n R 
^.0 ZIZ u.o 


0.91 


± 


0.27 




IVl \_y Vjr Zi 


1 q + n 4 

1 . y ^ U .4: 


Z( . y ziz u.o 


1 1 9 R 

1 . 1 ZIZ z . 


2.67 


± 


7.01 




Art 1 90 


o.z in u.o 


0.4 ziz U.O 


_|_ n 7 
0.0 ziz u. / 


0.95 


± 


0.12 




INVJV^ UU^O 


c 7 _i_ n 9 
0. 1 m u.z 


u.u ZIZ u.o 


O.U ZIZ u.o 


1.07 


± 


0.08 




/A.1J\. 0\J^ 


10. z HZ u.o 


1 Q n _|_ n 7 

10. U ZIZ U. 1 


99 f; _l_ 1 c 

ZZ.U ZC 1 . 


0.80 


± 


0.06 




NPP 7zlRQ 




O.U ZIZ ±.u 


Q _|_ 1 

U.y ZIZ 1.0 


1.26 


± 


0.27 




TRA'^ 1S'?9^-^Q9fi 
irL/^kj loozo oyzu 


1 K 7 -1- n 

10. 1 HI u.o 


1 c Q _|_ n 
i.o.u ZIZ u.o 


1 fi -I- 1 1 

lU.U ZC 1 . 1 


1.00 


± 


0.08 




IN vjrw ( 00 


1 9 4 + 
1Z.4 m u.o 


±u. 1 ZIZ u.u 


lu.z ZIZ u.o 


1.63 


± 


0.08 




IVlAll / uu 


i^.u HI u.o 


±0.0 ZIZ 1.1 


^A^ -1-99 

1^4:. 1 ZIZ Z . Z 


1.09 


± 


0.21 




iNVj^ ^tjyo 


y.z ziz u.o 


±u.z ZIZ u.o 


Q _|_ n Q 
y.u ZIZ u.y 


1.05 


± 


0.13 




NPP "1^1 fi 


u.u ziz u.z 


o.z ZIZ u.o 


0.0 ZIZ u.o 


1.42 


± 


0.14 




A/rpn_fi '^0 1 


1 7 q _l_ A 4 
1 / .y m u.^ 


91 7 _|_ n c 
zi . 1 ZIZ u.u 


1 c q _|_ A 7 

10.0 ZIZ U. I 


1.41 


± 


0.08 




NPP "^997 


1 9 G -1- fl 9 

IZ.U HI U.Z 


1 4 R -1- n 4 

X'l.O ZIZ u.^ 


11 0-1-0 4 
11. U ZC u.^ 


1.34 


± 


0.06 




NPP zlOt^l 


Q9 n + 1 9 

OZ .U HI 1 .Z 


4:U.Z ZIZ Z.U 


9f: c _l_ 9 f: 

zu.o ZC Z.U 


1.52 


± 


0.18 


fi H 

U Q 


■^P 1 90 


O.U HI U.O 


O.z ZIZ u.u 


0. / ZC u.u 


0.92 


± 


0.15 




P cLll dil a 


y zn 1 . z 


0.0 ZIZ z . z 


1 n 1 -1- 9 Q 
lu. 1 ZC z . y 


0.88 


± 


0.26 




TP /l'^9Qa 
iv^ ^ozya 


1 .0 m u.o 


7 7 _i_ n 
i . ( ZIZ u.u 


u. / ZC u.u 


1.15 


± 


0.14 




A/rii-n 97Q 
ivijvii z / y 


19 4-1-1 e. 

IZ.I HI l.U 


1 yl 1 -U Q 1 
i^t. J. ZIZ 0. -L 


« 9 _|_ /l Q 
u.z ZC ^.0 


2.27 


± 


1.66 




NPP ^^i4S 


1 9 « -L fi 

IZ.U HI U.U 


1 C _L 1 1 
±0.0 ZIZ ± . ±. 


1 9 _L 1 9 

IZ.U ZIZ 1 . z 


1.07 


± 


0.14 




ArV ^>f!zl 


Ql « -1- 1 9 

1 .U HI 1 .Z 


00 1 _|_ 1 c 

00.1 ZIZ 1.0 


on 9 ^ /I n 

OU.Z ZC ■4:.U 


1.10 


± 


0.15 




NPP 74fiQ 
iNv^^ ( ^uy 


1 9 Q -L n 4 

IZ.O HI U.^ 


10. U ZIZ U. ( 


11 ^ n Q 
11.0 ZC u.o 


1.10 


± 


0.10 




NPP '^78'^ 


1 9 a -1- n 

XZ.O HI U.O 


X'i.'i ZIZ ±.u 


1 9 1 _l_ 1 n 

IZ.l ZC l.U 


1.18 


± 


0.26 




A/TPP— fi '?0 1 


99 4 + 1^ 

ZZ.^ HI U.O 


97 1 -U n R 
z / . 1 ZIZ u.o 


1 Q n _|_ 1 n 
ly.u ZC l.U 


1.43 


± 


0.08 




NPP 41 1^1 


11 7-1-0 4 

1 1 . J HI U.^ 


10. Z ZIZ l.U 


11 ^ n 7 

1 1 .U ZC U. i 


1.13 


± 


0.11 




NPP "^997 


97 4 -L 1 Q 

Z f HI X .0 


00 1 -U 9 Q 
00. 1 ZIZ z. y 


9/1 9 _L 9 c 

Z1.Z ZC Z.U 


1.37 


± 


0.15 






ou.o HI u.y 


to.i zjz 1.0 


01 q _|_ 9 n 

01.0 ZC Z.U 


1.38 


± 


0.10 


'\f\ H 

OU Q 




97 c -L fi q 

Z / .0 HI U.O 


90 _|_ fi c 
zy.o ZIZ u.u 


9c Q _l_ n c 
ZO.O ZC u.u 


1.16 


± 


0.04 






c n + 9 

O.U m u.z 


R 4 -1- n 4 

0.4: ZIZ U.4: 


R 4 -1- ^ 
0.^ ZC u.o 


1.00 


± 


0.08 




J? cLll clil ij 


1 9 c, -L 1 9 

IZ .0 HI 1 .Z 


1 n -1- 9 9 

1 .U ZIZ z .z 


11 Q _|_ 9 7 
1 1 .y ZC z . / 


1.26 


± 


0.34 






11 7 -L n 

XX. 1 HI U.O 


iz. 1 ZIZ u.o 


1 n 7 -1- n p; 
lu. / ZC u.o 


1.18 


± 


0.07 




NPP '=i'=i4?^ 


9t; fi -L n t; 

ZO.U HI U.O 


9f; 9 _L fi Q 
zu.z ZIZ u.y 


9^ 4 _l_ 1 1 

ZtJ.rt 1 1.1 


1.03 


± 


0.06 






qq 9 -L 1 1 

00. Z HI J L 


00.0 ZIZ 1.^ 


qf? fl -1- *^ ^ 
ou.u ZC 0.0 


0.92 


± 


0.10 




NPP 7Af\Gk 


1 « n + 4 

lU.U HI U.^ 


1 ( .y ZIZ u.u 


1 4 ^ n R 

14:. U ZC u.o 


1.23 


± 


0.08 




IMVjtV^ OOlU 


1 q 7 -L n 4 

1 y . / HI u.^ 


9f5 c: _|_ n R 
zu.o ZIZ u.o 


1 A ^ n 7 

lU.'l ZC U. ( 


1.61 


± 


0.09 




IVlv^Vjr U-OU-l(J 


9c; c -L n 1^ 

ZO.O HI u.o 


ou.o ZIZ u.y 


99 7 _1_ 1 1 

ZZ. t ZC 1 . 1 


1.35 


± 


0.08 




NPP 41 ^1 


94 c; -L n 4 

Z^.O HI U.^ 


94 Q _l_ 1 n 
Z'i.y ZIZ l.U 


9/1 7 _i_ n 

Z'it. / ZC U.U 


1.01 


± 


0.05 




NPP "^997 


c;9 R -L Q 
oz.u HI U.y 


f^n _|_ 1 Q 
uu.o ZIZ 1.0 


/If? q _i_ 1 7 

^U.O ZC 1 . ( 


1.31 


± 


0.06 




NPP 4(1^^1 


41 q -L 1 n 

4: 1 .0 HI 1 .U 


^U.o ZIZ 1.0 


q^; 7 ^ 9 1 

OU. ( ZC Z. 1 


1.28 


± 


0.08 




pp ns^o/i-i-Vfii 

r\j UoU4-t-/Dl 


i^.Z Hi i.U 


1 9 -1- 1 (\ 

LO.Z Zt l.U 


1 n 9 -1- /I 

lU.Z ZC 0.4 


1.28 


± 


0.45 




3C 390.3 


22.0 ± 0.4 


25.6 ± 0.7 


20.3 ± 0.8 


1.26 


± 


0.06 




3C 120 


16.9 ± 0.4 


18.4 ± 0.7 


16.0 ± 0.8 


1.15 


± 


0.07 




Fairall 9 


22.7 ± 0.7 


24.6 ± 1.3 


20.8 ± 1.6 


1.18 


± 


0.11 




Ark 120 


20.9 ± 0.4 


22.9 ± 0.7 


18.8 ± 0.9 


1.22 


± 


0.07 




IC 4329a 


15.7 ± 0.2 


17.0 ± 0.4 


14.3 ± 0.4 


1.19 


± 


0.05 




Mkn 110 


53.2 ± 1.0 


57.4 ± 1.8 


49.9 ± 2.2 


1.15 


± 


0.06 




Mkn 590 


30.3 ± 0.8 


32.5 ± 1.2 


28.7 ± 1.5 


1.13 


± 


0.07 




NGC 5548 


27.2 ± 0.4 


28.9 ± 0.8 


26.3 ± 0.9 


1.10 


± 


0.05 




Ark 564 


27.7 ± 0.9 


29.1 ± 1.1 


29.1 ± 2.8 


1.00 


± 


0.10 




Mkn 335 


32.0 ± 1.1 


34.8 ± 1.5 


28.5 ± 2.7 


1.22 


± 


0.13 
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Table 3— Continued 


Time 


Source 


2-12 keV 


2-4 keV F^ar (%) 7 


12 kcV F,ar (%) 






Scale 


Name 


P D/if \ A) ) 


(F f ] 


(F ^. ^) 




r,so/t 




r-.hor-d 




Mkn 79 


26.0 ± 1.2 


27 A ± 2.3 


27.7 ± 2.6 


0.99 


± 0.12 




NGC 7469 


24.1 ± 0.6 


25.7 ± 1.1 


23.5 ± 1.3 


1.09 


± 0.08 




NGC 3783 


20.9 ± 0.2 


24.0 ± 0.4 


18.7 ± 0.5 


1.28 


± 0.04 




NGC 3516 


29.4 ± 0.3 


38.4 ± 0.6 


25.0 ± 0.6 


1.54 


± 0.04 




MCG-6-30-15 


25.7 ± 0.4 


30.4 ± 0.7 


22.2 ± 0.9 


1.37 


± 0.06 




NGC 4151 


32.3 ± 0.1 


37.7 ± 0.3 


30.0 ± 0.2 


1.26 


± 0.01 




NGC 3227 


34.9 ± 0.6 


41.0 ± 1.3 


30.4 ± 1.2 


1.35 


± 0.07 




NGC 4051 


61.7 ± 1.0 


76.0 ± 1.7 


51.8 ± 2.2 


1.47 


± 0.07 


izyo u 


Qr^ 1 on 


io.o =t U.O 


ZO.4 ± U.D 


i / . / =t U.D 




-I- n 
=t u.uo 




Fairall 9 


38.2 ± 0.8 


40.7 ± 1.2 


37.2 ± 1.5 


1.09 


± 0.06 




NGC 5548 


31.6 ± 0.4 


34.0 ± 0.7 


30.0 ± 0.8 


1.13 


± 0.04 




Ark 564 


37.1 ± 1.0 


42.5 ± 1.3 


32.4 ± 3.0 


1.31 


± 0.13 




NGC 3783 


20.6 ± 0.3 


22.7 ± 0.5 


17.6 ± 0.6 


1.29 


± 0.05 




NGC 3516 


32.0 ± 0.3 


40.0 ± 0.6 


25.4 ± 0.6 


1.57 


± 0.04 




MCG 6-30-15 


24.5 ± 0.3 


29.1 ± 0.6 


19.9 ± 0.8 


1.46 


± 0.06 




NGC 3227 


56.3 ± 0.6 


64.0 ± 1.2 


48.5 ± 1.1 


1.32 


± 0.04 




NGC 4051 


61.7 ± 0.8 


75.9 ± 1.4 


51.1 ± 1.7 


1.49 


± 0.06 



Note. — Targets are ranked by 2-12 keV luminosity, given in Table 1. 1 d F^^r for Mkn 279 is undefined 
at all bands and is not listed. 



Table 4 

Derived Variability Parameters for 



THE ASCA LIGHT CURVES 


Source 


2 lOkcV 


Nam(> 


.F,v„. m 


PG 0804-^761 


3.2 ± 1.9 


Mkn 509 


4.0 ± 0.7 


3C 120 


1.5 ± 0.9 


MCG-2-58-22 


3.1 ± 1.0 


Fairall 9 


3.4 ± 2.3 


Ark 120 


2.7 ± 1.1 


IC 4329a 


5.6 ± 0.7 


NGC 5548 


6.5 ± 0.9 


Ark 564 


17.0 ± 1.2 


NGC 7469 


2.2 ± 1.5 


IRAS 18325-5926 


14.0 ± 1.9 


NGC 3783 


6.1 ± 0.8 


Mkn 766 


17.6 ± 1.6 


NGC 4593 


14.1 ± 0.9 


NGC 3516 


10.3 ± 0.9 


MCG-6-30-15 


20.0 ± 1.2 


NGC 4151 


6.3 ± 0.4 


NGC 3227 


14.8 ± 0.9 


NGC 4051 


29.6 ± 2.0 
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Table 5 

Summary of Broadband Variability Correlations 



Figure 
Number 


X-axis 


y-axis 


Time 

Scale (d) 


No. 
Pts. 


r 


Pr 


Slope 


Xr 


2a 


1^2-12 






19 


—0.754 


1 HQ 1 n — 4 
i.yoXiU 


n A A _l_n mi 
— 0.344±U.011 


30.7 








1 \HX 1 a) 


17 


—0.690 


2.16X10 


— 0.316±U.024 


43.6 








6 


12 


-0.734 


6.54x10-3 


-0.223+0.012 


14.3 








36 


12 


-0.531 


7.54x10-^ 


-0.176+0.017 


15.4 








216 


19 


-0.642 


3 07x10-3 


-0 138+0 006 


6.4 








1296 


9 


-0.533 


0.139 


-0.135+0.014 


8.8 


2b 


Mbh 




1 (ASCA) 


18 


-0.670 


2.36x10-3 


-0.321+0.011 


41.7 








1 (RXTE) 


15 


-0.839 


9.14x10-^ 


-0.368+0.015 


18.7 








6 


12 


-0.441 


0.152 


-0.259+0.024 


23.1 








36 


12 


-0.098 


0.762 


-0.077+0.022 


30.1 








216 


19 


-0.396 


9.20x10-2 


-0.113+0.010 


10.4 








1296 


9 


-0.133 


0.732 


-0.046+0.026 


18.1 


3 


^var,3oft 


^var,hard 


1 


16 


+0.950 


1.83x10-® 


+0.868+0.030 


9.4, 50.8 




6 


12 


+0.867 


2.60x10-'' 


+0.800+0.021 


1.0, 26.8 








36 


12 


+0.923 


1.86x10-^ 


+0.960+0.017 


12.2, 21.9 








216 


19 


+0.925 


1.54x10-* 


+0.927+0.012 


9.4, 22.4 








1296 


9 


+0.967 


2.16x10-5 


+0.981+0.021 


17.8, 41.1 


4 


L2-I2 


^var,soft 


1 


16 


-0.585 


1.70x10-2 


-0.161+0.035 


6.6 




^var,hard 




















6 


12 


-0.571 


5.25x10-2 


-0.081+0.027 


1.4 








36 


12 


-0.476 


0.118 


-0.071+0.010 


6.1 








216 


19 


-0.358 


0.132 


-0.091+0.006 


6.5 








1296 


9 


-0.583 


9.92x10-2 


-0.103+0.012 


9.5 




Mbh 


^var,soft 


1 


14 


-0.033 


0.911 


+0.147+0.057 


10.7 




^var,hard 




















6 


12 


-0.445 


0.147 


-0.069+0.0.36 


1.6 








36 


12 


-0.084 


0.795 


-0.012+0.016 


6.4 








216 


19 


-0.074 


0.764 


-0.037+0.012 


6.7 








1296 


9 


-0.567 


0.112 


-0.126+0.014 


8.4 


6 


1/2-12 






23 


-0.700 


1.99x10-" 


-0.080+0.002 


310.26 


6 


Mbh 







23 


-0.690 


2.66x10-" 


-0.094+0.003 


277.62 



Note. — Col. (6) is the Spearman rank correlation coefficient. Col. (7) is the probability of obtaining that correlation 
coefficient by chance. The best-fitting slopes, listed in Col. (8), were obtained using the regression method of Akritas & 
Bershady (1996); their bisector regression was used for the correlations in Figure 3. The reduced chi-squaxed values Xr 
listed in Col. (9). The number of degrees of freedom used in calculating Xr is the number of data points (Col. [5]) minus 2 for 
all data sets. For all six data sets in Figures 2a and 2b, Xr calculated assuming a uniform error bar of 0.047 and 0.044 in 
the log, respectively (one-sided; the average error for all data points in the figure) . The first value of Xr fisted for the Figure 3 
correlations was calculated with F^ar,soft ^■s the independent variable and F^ar,hard as the dependent variable and vice-versa 
for the second value. For the 1 d '■ — - — A^bh correlation, the positive slope is caused by an outlier point, MCG-2-58-22. 

^ Tjar,hard 

Removal of this points yields r=-0.291, Pt.=0.334, a slope of -0.106+0.024, and x?=ll-0. Ark 564 and PC 0804+761 have 
been excluded from the second row of Figure 6 correlations. 
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Table 6 



Summary of color-flux diagrams 



Source 
Name 


MJD 
Range 


(HR) 









IC 4329a 


52737- 


-53035 


1 


.10 


1, 


.178±0, 


.008 


NGC 4151 


51678- 


-51964 


o 

z 


.Do 


1, 


.075±0, 


.003 


NGC 3783 (2) 


519 m- 


C 00*70 

-5zz ( z 


-1 

1 


.16 


1 


.207±0, 


.032 


NGC 3783 (1) 


516 


-51975 


-1 

1 


1 n 
.19 


1 


.188±0, 


.009 


3C 120 


oz3oo- 


c oc o c 


1 


.08 


1, 


.042±0, 


.005 


MCG-6-30-15 (3) 


52276- 


-52554 





.93 


1, 


.308±0, 


.010 


MCG-6-30-15 (1) 


51676- 


-51974 


0, 


.97 


1 


.224±0, 


.015 


MCG 6-30-15 (2) 


51976- 


-52274 





.98 


1, 


.224±0, 


.006 


NGC 5548 (2) 


51977- 


-52271 


1 


.06 


1, 


.032±0, 


.003 


NGC 3516 (2) 


50838- 


-51119 


1 


.43 


1, 


.318±0, 


.017 


NGC 3227 (2) 


52300- 


-52599 


1 


.31 


1, 


.269±0, 


.019 


NGC 5448 (1) 


51678 


51968 


1 


.06 


1, 


.024±0, 


.006 


NGC 3516 (1) 


50523- 


-50799 


1 


.40 


1, 


.388±0, 


.004 


NGC 3227 (1) 


52000- 


-52298 


1 


.60 


1, 


.238±0, 


.011 


3C 390.3 


51261- 


-51555 


1 


.04 


1, 


.059±0, 


.004 


Ark 120 


50868- 


-51167 





.97 


1, 


.075±0, 


.018 


NGC 4051 (3) 


52276- 


-52574 


1 


.00 


1, 


.309±0, 


.010 


NGC 4051 (1) 


51676- 


-51974 


1 


.03 


1, 


.303±0, 


.005 


NGC 5548 (3) 


52276- 


-52574 


1 


.09 


1, 


.016±0, 


.022 


Fairall 9 (2) 


51978- 


-52268 





.97 


1, 


.116±0, 


.009 


Fairall 9 (3) 


52277- 


-52571 


1 


.00 


1, 


.162±0, 


.043 


NGC 4051 (2) 


51976- 


-52274 


1 


.03 


1, 


.241±0, 


.010 


Ark 564 (1) 


51678- 


-51973 





.43 


0, 


.921±0, 


.008 


Ark 564 (2) 


51977- 


-52271 





.47 


0, 


.950±0, 


.056 


Ark 564 (3) 


52276- 


-52574 





.50 


0, 


.861±0, 


.044 


Fairall 9 (1) 


51679- 


-51974 


1 


.03 


1, 


.057±0, 


.010 


PG 0804+761 


51676- 


-51971 





.71 


0. 


.813±0, 


.026 



Note. - Targets are ranked by mean geometric mean count 
rate. Col. (3) is the mean 7-12 keV / 2-4 keV count rate hardness 
ratio (HR) value for light curve, obtained from the binned data in 
Figure 6. Col. (4) lists <p, which quantifies the decrease in HR for 
every doubling in flux, as described in the text. 
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Table 7 

Summary of Count Rate Correlations 



Figure 






Time 


No. 






Number 


X-axis 


y-axis 


Scale 


Pts. 


r 


Pr 



8 Count rate 


Fvar 


1 d {ASCA) (all srcs.) 


48 


-0.240 


0.100 






1 d (RXTE) (all srcs.) 


84 


-0.192 


0.080 






6 d (all srcs.) 


67 


-0.516 


7.99x10"*' 






36 d (all srcs.) 


19 


-0.577 


9.67x10-3 






216 d (all srcs.) 


78 


-0.415 


1.60x10-* 






1296 d (all srcs.) 


12 


-0.816 


1.21x10-3 






1 d (ASCA) (BL only) 


31 


-0.202 


0.276 






1 d (RXTE) (BL only) 


56 


-0.308 


0.021 






6 d (BL only) 


39 


-0.479 


2.06x10-3 






36 d (BL only) 


14 


-0.481 


0.083 






216 d (BL only) 


53 


-0.202 


0.147 






1296 d (BL only) 


7 


-0.786 


0.036 






6 d (BL only, NGC 3227 cxcl.) 


29 


-0.217 


0.258 






36 d (BL only, NGC 3227 excl.) 


12 


-0.245 


0.443 


- ^2-12 


Count rate 


1 d (ASCA) 


21 


-0.161 


0.487 






1 d (RXTE) 


17 


+0.158 


0.544 






6 d 


12 


+0.070 


0.829 






36 d 


12 


+0.042 


0.897 






216 d 


19 


-0.229 


0.346 






1296 d 


9 


+0.218 


0.574 


Mbh 


Count rate 


1 d (ASCA) 


19 


-0.011 


0.966 






1 d (RXTE) 


15 


+0.378 


0.164 






6 d 


12 


-0.021 


0.948 






36 d 


12 


-0.035 


0.914 






216 d 


19 


-0.063 


0.797 






1296 d 


9 


+0.167 


0.667 


9 Geom. mean count rate 


^var,soft 
^var,hard 


1 d 


80 


+0.499 


2.14x10-^ 






6 d 


67 


+0.375 


1.50x10-3 






36 d 


19 


-0.077 


0.753 






216 d 


78 


+0.472 


1.30x10-5 






1296 d 


12 


+0.077 


0.812 


10 Geom. mean count rate 


<t> 




27 
23 


+0.405 
+0.056 


3.59x10-2 
0.799 



Note. — Col. (6) is the Spearman rank correlation coefficient. Col. (7) is the probability of obtaining that correlation 
coefficient by chance. For the Figure 10 correlations, Ark 564 and PC 0804+761 have been excluded from the second row. 
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Fig. 1. — RXTE 2-12 keV total light curves, ranked by 2—12 keV luminosity, before clipping and resampling. The black, red, green, blue 
and orange bars denote the extent of the 1 d, 6 d, 36 d, 216 d, and 1296 d RXTE light curves, respectively, before subsampling to a common 
sampling rate. Error bars are la. 
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Figure 1, cont'd. 
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Fig. 1. — Figure 1, cont'd. 
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Fig. 2. — Fyar plotted against 2-10 keV luminosity (left) and black hole mass estimate A/bh (right) for all time scales. 1 d ASCA, 1 d 
RXTE, 6 d 36 d, 216 d, and 1296 d time scale data points are denoted by purple open circles, black filled circles, red open squares, green 
filled squares, blue open triangles, and orange filled triangles, respectively. The best-fit lines for each time scale are the same color. 




4 F (%) 



var 

Fig. 3. — 7-12 keV Fyar plotted against 2-4 keV F^ar- Data points and best-fit lines are denoted the same as in Figure 2. A source with 
equally strong variability in the two bands would lie on the dashed line, but the vast majority of the light curves exhibit stronger variability 
in the softer band. 
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Fig. 4. — The ratio of 2-4 keV F^ar / 7-12 keV F^ar plotted against source luminosity. Data points and best-fit lines are denoted the same 
as in Figure 2. There is tentative evidence for relatively less luminous sources to display increasingly stronger variability in the soft band 
compared to the hard band. 



26 



Seyfert 1 X-ray Variability 




Fig. 5. — Plots of the 7-12 keV/2-4 keV hardness ratio (HR) against geometric mean count rate for the 27 hght curves with adequate data. 
Gray error bars represent the unbinned data; black points represent the binned data. The solid lines are the best linear fits to the binned 
data. 
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Fig. 5. — Continued. 
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Fig. 6. — The parameter (j>, which quantifies the decrease in 7-12 keV/2-4 keV hardness ratio (HR) for every doubUng in fiux, plotted 
against 2—12 keV luminosity and black hole mass estimate Mbh- Ark 564 and PG 0804+761 are excluded. There is a tendency for relatively 
less luminous or less massive sources to display more overall spectral variability. 
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Fig. 7. — Predicted and measured ratios of Fvar for six time scale combinations plotted against black hole mass estimate A/bh (left) and 
2-12 keV luminosity (right). The solid lines and dotted lines denote the model Fyar ratios, derived from linearly scaling the PSD break 
frequency with black hole mass using the best-fit relation Tj,=A/bh/10®'^ Mq. The model ratios for the Fyar ratio — L2-12 plot were derived 
from linearly scaling PSD break frequency with X-ray luminosity using the best-fit relation Ti,=L2— 12 ^rg s~^. The black, red, 

green, and blue solid lines mark the model Fvar ratios -Ri/g, R1/3Q, R1/216 '^.nd fli/12961 respectively, where e.g., Ri/g denotes the ratio 
Fyar(^ d) / Fvar{5 d). The ratios -R6/2161 ^6/1296 3.re marked by cyan and orange dotted lines, respectively. Note that, as one scales a PSD 
towards higher temporal frequencies, Fvar measured over a fixed time scale (fixed sampling window in the frequency domain) will increase, 
but at a faster rate for relatively shorter sampling time scales. The model ratios hence follow iJi/g > ^1/36 > ^1/216 > ^1/1296 {^^^ black, 
red, green, and blue curves, respectively), -R5/36 > ^3/216 > ^6/1296 {''^ot plotted, cyan, and orange, respectively), etc. Note also that the 
Fvar ratio will tend to fiatten off at high frequencies when both sampling time scales are much shorter than the PSD break time scale; i.e., this 
is why i?i/g flattens off at high mass/ low frequencies. The measured ratios of Ri/g, ^2i/36i ^1/216> ^1/I296i ^6/2161 ^6/l296i ^.re denoted 
by black filled circles, red open triangles, green filled squares, blue crosses, cyan open circles and orange filled triangles, respectively. 
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Flo. 8. — Logarithm of 2—12 kcV F^ar (2-10 kcV for ASCA) plotted against the logarithm of the count rate. Filled circles denote broad-line 
Seyferts. Open circles denote narrow-line/ soft spectrum Seyferts, which tend to be more variable and less luminous in the RXTE bandpass. 
Filled triangles in the 6 d and 36 d plots axe NGC 3227; excluding those data points, no trends with count rate axe evident for the broad-line 
sources overall. 
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Fig. 9. — F^ar,soft/ Fvar,hard plotted against the logarithm of the count rate. No trends with count rate are evident. 
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Fig. 10. — Plot of the parameter (f) against average geometric mean count rate. Open circles denote Ark 564 and PG 0804+761; their 
low values of </> may be an artifact of the low count rates and systematic errors in the PCA background model. For the other sources, </> is 
independent of count rate. 



